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INTRODUCTION 


This  report  covers  progress  made  in  the  period  1  April,  1984  to 
31  March,  1985  for  the  eleven  research  units  funded  under  the  Joint  Services 
Electronics  Progr^un  at  Harvard  University.  It  is  broken  down  into  four 
major  divisions  of  electronic  research — solid  state  electronics,  quantum 
■electronics,  information  electronics  control  and  optimization,  and 
electromagnetic  phenomena.  It  also  includes  a  report  of  Significant 
Accomplishments  which  contains  noteworthy  results.  Ihese  are: 

X-Ray  Reflectivity  Measurement  of  the  Structure  of  the  Surface  of  Water 
by  Research  Unit  3,  Observation  of  Chaos  in  Josephseon  Junctions  Driven 
at  Far- Infrared  Frequencies  by  Research  Unit  4,  and  The  Application  of 
Electromagnetic  Surface  Waves  to  Geophysical  Prospecting  from  the  Surface 
of  the  Earth  by  Research  Unit  11.  This  significant  acconplishment  report 
is  not  duplicated  in  other  sections  of  the  Annual  Report  where  other 
activities  funded  under  the  same  research  unit  are  described. 
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1.1  Electronic  Structure  of  HqMnSe  and  Related  Semiconducting  Alloys 


H.  Ehrenreich,  K. C-  Hass,  amd  R.  Lempert,  Ctontract  N00014-84-K-0465; 
Research  Unit  I. 


!Ihis  group  has  been  continuing  work,  '  reported  in  previous  proposals 
and  euinual  reports,  concerning  the  application  of  a  semi-empirical  tight- 
binding  approach  to  the  calculation  of  the  electronic  structure  cund  other 
properties  of  2-6  alloys.  During  the  past  year  this  effort  has  emphasized 
alloys  containing  Mn  which  have  interesting  magnetic  properties  of  possible 
practical  importeuice.  The  peureuneters  required  for  in^lementation  of  the 
empirical  tight-binding  scheme  are  obtained  from  e:q>erimental  data  and 
local  spin  density  ASW  calculations  performed  in  collaboration  with  A.E. 
Carlsson  (University  of  Washington,  St.  Louis).  Ihese  calculations  pertain 
to  a  hypothetical  tetradxedral  antiferromagnetic  NnTe  and  MnSe.  The  same 
type  binding  techniques  are  also  being  utilized  within  the  coherent  potential 
approximation  for  a  calculation  of  the  electronic  structure  of  a  pseudo- 
quartemary  alloy  consisting  of  Cd,  Mnf,  Mn-f,  eund  Te.  The  results  will 


permit  examination  of  the  effects  of  disorder  on  both  the  magnetic  and 
electronic  properties  of  the  alloy.  Correlation  effects  eu:e  included  within 
the  Hiihdeurd  model  amd  treated  within  the  Hcurtree-Fock  approximation.  The 
various  approaches  used  here  together  with  the  veuriety  of  semiconducting 
dilute  magnetic  alloys  being  examined,  should  provide  a  reliable  overview 
of  the  electronic  structure  of  these  materials. 
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Hgi.jjMtijjTe  and  Hg^.jfOdxTe , "  Proceedings  of  the  MCT  Workshop,  J.  Vaa. 
Sot.  Teohnol.  Al,  1678  (.1983). 

2.  K. C.  Heiss,  B.  Velicky,  and  H.  Ehrenreich,  "Simplification  of  Green's 
Function  Calculations  Through  Analytic  Continuation,"  Phya.  Rbv.  B29, 
3697  (1984) . 


1.2  Magnetic  Properties  of  Diluted  Magnetic  Semiconductors.  H.  Ehrenreich, 

R.  Lenpert,  and  K.C.  Hass,  Contract  N00014-84-K-0465;  Research  Unit  1. 

Because  of  the  unusual  magnetic  properties  of  dilute  magnetic  semi¬ 
conductors,  there  heis  been  a  resurgence  of  interest  in  various  fundamental 
aspects  of  the  theory  of  their  magnetism.  These  pertain  to  the  exchange 
interaction  between  the  Mn  moments  and  conduction  bemd  electrons  or  valence 
bemd  holes  eis  well  as  the  interaction  among  Mn  moments  themselves.  The 
exchange  associated  with  the  electrical  carriers  is  reponsible  for  the  pro¬ 
nounced  and  possibly  useful  magneto-optical  effects.  That  associated  with 
the  Mn-Mn  interactions  is  inportant  in  the  understanding  of  cluster  anti¬ 
ferromagnetism  and  the  spin-glass  behavior  that  is  observed  at  low  tempera¬ 
tures.  We  are  currently  examining  the  hierarchy  of  Mn-Mn  exchange  mechanisms 
in  2-6  alloys  such  as  those  mentioned  above  using  a  generalized  Anderson 
Hamiltonian  and  the  electronic  structure  already  described.  The  hybridization 
strength  is  extracted  from  experimental  sp-d  exchange  constant*  using  an 
analytic  expression  provided  by  the  Schrieffer-Wolff  transformation.  All 
contributions  to  the  Mn-Mn  exchange  for  nearest  euid  second  nearest  neighbors 
aure  calculated  perturbative ly  to  lowest  non-vanishing  order.  Superexchemge 
is  found  to  dominate  in  agreement  with  the  observed  amtiferromagnetic 


coupling.  Other  contributions  including  the  Bloenibergen- Howland  mechanism 
and  direct  exchange  are  more  than  an  order  of  magnitude  smaller.  Hie  agree¬ 
ment  between  the  calculated  and  experimental  values  of  the  nearest  and  next 
nearest  neighbor  exchange  constants  is  very  satisfactory.  Similar  calculations 
for  other  dilute  magnetic  semiconductors  are  in  progress.  They  suggest  that 
the  dominance  of  superexchange  is  a  general  feature  of  this  class  of  alloys. 

The  chemical  trends  associated  with  changes  in  the  non-roagnetic  constituents 
are  also  under  examination  by  using  these  results  and  siirplified,  physically 
transpeurent  perturbation  expressions.^ 
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1.3  Semiconducting  Pseudobinary  Alloys.  H.  Ehrenreich,  K. C-  Hass,  and  R. 

Len^ert,  Contract  N00014-84-K-0465;  Research  Unit  1. 

The  recent  experimental  discovery  of  a  bimodal  distribution  of  nearest 

neighbor  bond  lengths  in  (In,  Ga)As  and  Zn(Se,Te)  has  raised  new  questions 

concerning  the  role  of  structural  disorder  in  semiconducting  alloys.  IXiring 

the  past  year  we  have  presented  the  first  self-consistent  effective  medium 

theory  of  the  electronic  structure  of  pseudobinary  zincblende  alloys  and 

applied  it  to  both  alloy  systems.^  These  calculations  rely  crucially  on  the 

molecular  coherent  potential  approximation,  which  can  be  readily  applied  to 

pseudobinary  alloys  euid  includes  the  effects  of  both  chemical  euid  structural 

disorder.  The  results  are  surprising:  besides  showing  that  the  effects  of 

structural  and  chemical  disorder  are  comparable  in  magnitude,  there  is 

evidence  that  both  disorder  effects,  when  considered  together,  can  either 

enhance  or  diminish  the  total  scattering  of  electrons  in  particular  energy 

regions.  The  mobility  of  conduction  band  electrons  is  actually  enhanced  in 

the  systems  under  consideration.  The  question  of  why  these  alloys,  in  view 

of  the  very  appreciable  associated  lattice  distortion  should  not  form  in  the 

2 

chalcopyrite  structure  remains  open  and  is  under  investigation. 
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1.4  Dis1ocation~Free  GaAs.  H.  Ehrenreidi,  Contract  N00014-84-K-0465;  Research 
unit  1. 

This  research,  performed  in  collaboration  with  J.P.  Hirth  (Chio  State 
University)  examines  the  question  of  why  the  addition  of  a  small  ^uaount  of  In 
to  as- grown  GaAs  crystals  should  result  in  2m  appreciable  reduction  in  disloca¬ 
tion  density.  The  model  proposed  to  explain  this  phenomenon  is  based  on  solid 
solution  heurdening,  familiar  in  metallurgy.  In  the  present  case  appreci2d}le 
differences  in  the  volvime  fraction  to  be  ckssociated  with  solute  emd  solvent 
atoms  respectively  cure  likely  to  produce  substamtial  reductions  in  dislocation 
density  of  the  host  crystal  without  deleterious  effects  on  the  electrical 
properties.  Ihis  feature  is  likely  to  be  relevemt  to  alloy  systems  other  than 
InGaAs.  Furthermore,  it  may  possibly  provide  some  insist  concerning  the 
stability  of  strained  superlattices.  Ihe  results  will  appear  shortly  in 
Applied  Phyaioa  Letters,^ 

Refe rente 

1.  H.  Ehrenreich  and  J.P.  Hirth,  Appl.  Phya.  Lett. ,  to  be  published  ;^ril  1, 
1985  issue. 


1.5  Charged  Dislocations  and  Jogs  in  HgCdTe.  H.  Ehrenreich,  Contract 
N00014-84-K-0465f  Research  Uhit  1. 

Jogs  in  dislocations  in  2-6  coitpounds  act  as  extrinsic  sources  for 
vac2mcies  emd  interstitials.  They  carry  effective  chfurges  and  will  therefore 
act  2m  donors  emd  acceptors.  As  a  result  they  cem  influence  the  electrical 
properties  of  these  conpounds  signific2mtly.  This  research,  performed  in 
coll2d>oration  with  J.P.  Hirth  (Chio  State  University)^  suggests  reasons  for 
the  fact  that  HgOdTe  alloys,  which  would  ordin2u:ily  be  p-type  at  low 


dislocation  densities  converts  to  n-type  after  cleavage  on  a  [110]  plane  or 
mechanical  polishing.  Ihe  insight  provided  is  of  considerable  practical 
inportemce  euid  probably  of  more  general  applicability. 


Reference 


1.  J.P.  Birth  euid  H.  Ehrenreich,  "Charged  Dislocations  and  Jogs  In  Hg^.^Od^Te 
and  Other  II- IV  Ootnpoiinds ,  ”  J.  Vac.  Sai,  Teeth,  ,  to  be  published  in 
Naorch/ April  1985  issue. 


1 . 6  Optical  Moment  Sum  Rules.  H.  Ehrenreich  and  K. C-  Hass ,  Contract 
N00014-84-K-0465;  Besearch  Unit  1. 

A  methodology  is  introduced  which  utilizes  relationships  between  moments 
of  ordered  and  disordered  systems  as  a  test  of  the  validity  of  elec¬ 

tronic  structure  and  chemical  bonding  models.  The  zeroth  moment,  related 
within  a  one-electron  approximation  to  an  average  dipole  matrix  element,  is 
particularly  useful.  We  have  made  a  detailed  conparison  of  optical  data  for 
a  v^u:iety  of  eunorphous  semiconductors  and  their  corresponding  crystals. 

Despite  limited  avedlable  data  the  present  results  strongly  suggest  that  the 
difference  between  average  dipole  matrix  elements  associated  with  the  amorphous 
and  crystalline  phase  respectively  is  typically  less  than  15%  for  a  given  semi¬ 
conductor.  The  difference  is  positive  for  tetrahedral  semiconductors,  zero  for 
tetrahedral  oxides,  auid  negative  for  moleculeir  semiconductors.  These  results 
have  been  used  to  discuss  the  validity  of  the  "independent  band  model. "  Their 
significance  has  been  explored  by  utilizing  a  real  space  picture  which  focuses 
on  the  effects  of  the  short  remge  environment.  Some  questions  concerning  the 
structure  of  amorphous  semiconductors  can  be  clarified  with  the  help  of 
analyses  such  as  that  developed  through  this  research.^ 

Reference 

1.  K. C.  Hass  and  H.  Ehrenreich,  "Electronic  Structure  Models,  Bonding,  and 
Optical  Moments  in  Amorphous  and  c:rystalline  Semiconductors,"  kmals  of 
Physios,  to  be  published  in  October  1965  issue. 


1.7  Photovol talcs.  H.  Ehrenreich,  Contract  N00014-84-K-0465;  Research  Unit  1. 

During  the  period  1978-1980  the  senior  investigator  was  Chairroeui  of  a 

Nhite  House  Panel  on  Photovoltaic  Energy  Sources  which  was  sponsored  by  the 

Americam  Physical  Society.  Since  then  he  has  served  on  a  number  of  government 

and  industrial  panels  that  monitor  the  evolution  of  this  technology.  Most 

recent  views  aure  expressed  in  a  paper  in  Solar  Celle^  auid  in  testimony  before 

the  Subcommittee  on  Energy  Development  auid  Applications  of  the  Committee  of 
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Science  and  Technology,  House  of  Representatives,  Maurch  5,  1985. 

References 
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1.8  Investigation  of  Amorphous  Hydrogenated  Si  and  Si-Ge  Alloys.  J-l.  Hanna, 
P.B.  Kirby,  D.  Leopold,  K.  Mackenzie,  Z.  Sun,  J.  Burnett,  J.  Eggert,  Y-H 
Li,  and  W.  Paul,  Contracts  N00014-84-K-0465,  NSF  DMR-81-08327  and  SERI 
Subcontract  XB-2-02144-1  of  DOE  Prime  Contract  DE- AC-02-83- 0110093;  Research 
Unit  2  (former  #3). 

During  this  period  we  con^leted  our  investigation  of  the  properties  of 
amorphous  hydrogenated  alloys  of  Si  auid  Ge  made  by  glow  discharge  decomposition 
of  appropriate  mixtures  of  SiH^  amd  GeH^.  These  alloys  have  smaller  forbidden 
energy  gaps  than  a-Si:H,  which  potentially  make  them  more  suitable  for  most 
applications.  In  practice,  they  are  usually  found  to  have  inferior  photo- 
electronic  properties,  whose  cause  is  still  under  investigation.  The  alloys 
investigated  were  made  under  conditions  of  different  substrate  ten^erature, 
r.f.  power,  and  gas  tenperature  auid  flow.  After  the  usual  chemical  and 
structural  chauracterization  was  done,  measurements  were  made  of  electrical 
conductivity,  optical  absorption,  photoconductivity  and  photoluminescence. 

These  measurements  were  then  correlated  in  order  to  give  a  self-consistent 
model  of  the  electronic  band  structure  auid  properties.  A  lengthy  account  of 
the  systematic  investigation  of  these  alloys  made  from  the  deconposition  of 
SiH^-GeH^  mixtures  hais  been  accepted  for  publication  in  the  Phyaioal  Beview. 


During  this  period  we  have  also  undertaken  to  produce  amorphous  alloys 
of  Sij^_j^Ge^  by  the  glow  dischuge  decouposition  of  mixtures  of  SiF^,  GeF^  and 
The  idea  of  vising  F  as  a  compensator  atom  for  these  alloys  was 
suggested  four  years  ago  by  the  Principal  Investigator,  in  order  to  improve 
the  microstructure  of  the  £d.loys  by  eliminating  weak  bonds.  Subsequently  a 
Japainese  group  discovered  that  material  made  in  this  fashion  had  superior 
photoconductive  properties.  We  have  carried  out  the  same  type  of  investigation 
for  alloys  made  from  these  materials  as  from  the  hydrides,  vaurying  the  gas 
pressures  and  flow  rates  and  the  applied  r.f.  power.  The  same  battery  of 
investigative  property  measurements  was  done.  Ihe  result  at  the  time  of 
writing  is  that  the  photoconductive  response  is  improved  by  an  order  of  magni¬ 
tude.  Other  properties  have  also  changed  to  some  extent.  The  interesting 
fundamental  aspect  of  this  work  is  that  it  appears  likely  that  the  disorder 
in  the  material  may  have  been  altered  by  the  difference  in  compensation  method. 
Thus  there  is  a  possibility  that  the  extent  of  the  bcuid  tails  has  been  altered 
and  that  a  different  kind  of  diamagnetic  defect  state  is  prevalent  in  this 
material  from  that  obtained  in  the  alloys  produced  from  the  hydrides.  The 
implications  for  the  structure  amd  properties  of  Si-Ge  alloys,  and  other  two- 
component  alloys  with  compensated  dangling  bonds,  aure  still  under  investigation. 


1.9  Staebler-Wronskj  Investigation  of  a-Si :H,  Doped  a-Si  :H,  a-Slj_^6e^:H 

Alloys  and  (P.B)  Doped  and  Compensated  a-Si :H.  D.J.  Leopold  and  W.  Paul, 
Contracts  NSF  DMBr 81-08327  and  N00014-84-K-0465;  Research  Ikiit  2  (former  #3) 

The  outstanding  problem  in  the  full  utilization  of  a-Si:H  amd  related 
amorphous  materials  for  photoelectronic  applications  is  the  decreaise  in  photo¬ 
conductive  response  as  a  result  of  prolonged  light  irradiation.  It  is  commonly 
believed  that  the  cause  is  the  production  of  a  defect  (dangling-bond)  state  in 
the  recombination  event  which  follows  light  Irradiation.  We  have  been  in¬ 
vestigating  the  production  of  such  defects  by  measuring  the  changes  in  sub¬ 
band  gap  absoaption  euid  electrical  transport  caused  by  measured  doses  of 
radiation.  The  types  of  sample  used  have  been  a-Si;H,  a-Si^_^Ge^:H  and  doped 
emd  compensated  a-Si:H.  The  excitation  spectral  dependence  of  the  changes  in 
transport  and  optical  absorption  have  been  measured. 


The  particular  thrxist  of  our  recent  measurement  heis  been  the  differences 
in  properties  between  singlydoped  a-Si:H  and  doped/compensated  material. 

Using  the  constemt  photociirrent  method  to  extend  direct  optical  absorption 
measurements  to  sub-band  gap  energies ,  we  have  shown  that  doping  a-Si :H  either 
n-type  with  phosphorus  or  p-type  with  boron  results  in  an  Increued  defect 
density  neeu:  midgap.  Oos^naated  a-Si:H  samples  show  a  large  reduction  in 
defect  ^  nsity  compared  with  that  in  both  singly-doped  and  undoped  material. 
Also,  the  steady  state  photoconductivity  is  found  to  decrease  when  boron-d^ed 
material  is  compensated  with  phosphorus.  Tnls  is  consistent  with  time-of- 
f light  measurements,  idiich  show  a  reduced  drift  mobility  for  elect -ons  and 
holes  in  coiq>en8ated  material,  implying  the  existence  of  new  defe^  _  states  in 
the  band  tails.  The  absorption  spectra  st^port  the  proposal  that  additional 
defect  states  are  created  in  the  band  tails  as  increasing  asiounts  of  dopant 
are  incorporated  in  caiQ>ensated  a-Si:H.  The  observed  band  tail  brw^denlng 
presumably  reflects  the  presence  of  these  additional  defect  states. 

Light  soalclng  experiments  performed  on  doped  a-Si:H  sanples  show  a  dark 
conductivity  decrease  which  is  ooo^letely  reversible  upon  annealing.  Also, 
changes  in  fhotoconductivity  recombination  from  a  bimolecular  to  a  monomolecular 
controlled  process  have  been  observed.  03a%>ensated  material  can  give  rise  to 
either  a  dark  conductivity  decrease  or  increase  aifter  eiqposure  to  light, 
depending  on  the  relative  amounts  of  phosphorus  and  boron  incorporated  in  the 
material.  We  have  found  activation  energy  shifts  concomitant  with  these  duk 
conductivity  changes  induced  during  light  soaking.  Presumably  the  Fermi  level 
position  moves  in  response  to  additional  defect  states  created  during  exposure. 
Direct  evidence  for  new  defect  states  in  the  gap  wm  inferred  from  a  low  energy 
absorption  Increeise  after  optical  ejqjosure. 


1.10  Hetastable  Semiconducting  Alloys  of  Gci.^Sny.  S.  Lee,  J.  Burnett,  and 
W.  Paul,  Oontract  N00014-84-K-0465;  Besearch  Unit  2  (former  #3). 


We  have  continued  to  make  amorphous  Ge-Si  alloys  by  r.f.  sputtering  in 
an  Ar  plasma  of  a  polycrystalline  Ge  target  to  thlch  pieces  of  polycrystalline 
6-Sn  are  attached.  The  composition  of  the  resulting  films  is  determined  by  an 
electron  microprobe,  crystallization  is  accomplished  by  melting  with  one  10  ns 
pulse  from  em  excimer  latser,  in  collaboration  with  personnel  of  Lincoln 
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Laboratory,  emd  thereafter  X-ray  diffraction,  Reunan  and  electroreflecteuice 
studies  are  made.  Since  our  last  report,  we  have  concentrated  on  understanding 
and  separating  out  artifacts  in  the  electroreflectivity  experiment,  caused  by 
possible  interferences  due  to  the  polyimide  insulating  film,  from  true  critical- 
point  electroreflectivity  features  in  the  crystallized  Ge-Sn  alloy. 

That  there  w2is  indeed  the  possibility  of  interference  effects  in  the 
electroreflectivity  spectrum  was  verified  by  measuring  the  electrically- 
unmodulated  reflectivity  spectrum  of  crystalline  Ge  on  which  the  usual  thick¬ 
ness  of  polyimide  had  been  deposited.  Maxima  and  minima  in  the  reflectivity 
were  found  at  the  same  energies  as  in  the  electroreflectivity  spectrum,  except 
for  the  well-established  true  critical-point  features  at  2.1  auid  2.3  eV.  Since 
the  method  of  deposition  of  the  polyimide  film  gives  a  very  uniform  film  of 
about  the  same  thickness  in  each  run,  it  seemed  reasoncd^le  that  at  least  some 
of  the  maxima  and  minima  in  the  electroreflectivity  spectra  of  c-Ge,  c-Si,  and 
c-Gej^__^Sn^  were  related  to  interferences  in  the  polyimide.  The  photon  energies 
of  the  observed  maxima  amd  minima  were  consistent  with  the  optical  thickness  of 
the  polyimide  film.  However,  since  the  eiectroreflectivity  spectrum  measured 
only  the  change  in  the  reflectivity  when  an  electric  field  is  applied,  it  is 
necessaury  that  the  field  change  the  optical  constants  of  the  semiconductor- 
polyimide  package.  It  is  presumed  that  the  most  likely  cause  of  this  is  a 
chamge  in  the  refractive  index  of  the  semiconducto  r,  but  since  a  study  of  this 
is  not  the  maun  objective  of  our  investigation,  we  have  tried  instead  to  avoid 
the  phenomenon.  It  should  be  remarked  in  passing,  however,  that  these  inter¬ 
ference  effects  in  the  eiectroreflectivity  spectrum  aure  not  observed  for  the 
amorphous  Ge^  Sn  layers,  while  they  appear  for  the  crystallized  Ge-_  Sn  . 

This  would  appear  to  indicate  a  difference  in  the  effects  of  am  electric  field 
on  the  refractive  indices  of  these  amorphous  and  crystallized  semiconductors. 

In  principle,  these  interference  effects  may  be  eliminated  by  roughening 
the  crystal  film  surface  before  depositing  the  polyimide.  h  comparison  wais 
made  of  the  eiectroreflectivity  spectra  of  two  codeposited  Ge  films,  one  of 
which  had  been  etched  with  CP4  or  H2O2 d) :HI’(2) ;H20(30)  before  deposition  of 
the  polyimide  layer.  This  film  showed  augmented  eiectroreflectivity  features 
at  the  energies  of  the  (2.1  eV)  and  +  (2.3)  tramsitions,  amd  no 

other  oscillations.  When  a  film  of  crystallized  Ge^  ^Sn^  ^  etched,  all  of 
the  features  observed  in  the  eiectroreflectivity  spectrum  of  am  unetched  film 
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(HASYLAB)  in  Hamburg,  Germany.  The  most  important  technical  development  that 
went  into  this  spectrometer  was  a  channel  cut  crystal  monochromator,  Ge(lll), 
in  which  the  beam  is  selected  by  two  Bragg  reflections.  The  monochromator  is 
oriented  such  that  the  first  reflection  bends  the  beam  downweurd  by  the  Bragg 
angle  (29)  of  27.2  degrees.  The  second  reflection  bends  the  beam  upweurd  by 
am  identical  angle  such  that  the  emitted  monochromatic  beam  is  paurallel  to  the 
incident  beam.  The  monochromator  and  associated  slits  aure  mounted  on  a  rigid 
bar  whose  position  and  orientation  are  cooputer  controlled.  In  this  way  a 
single  well-defined  downwaurd  incident  angle  is  selected  from  the  output  of  the 
rotating  amode  x-ray  source.  The  sanple  is  on  a  conputer  controlled  elevator 
that  is  continually  adjusted  so  that  the  monocnromatic  beam  strikes  the  center 
of  the  sample.  The  reflected  beam  is  detected  by  a  scintillation  detector-slit 
combination.  Using  this  spectrometer  it  was  possible  to  make  preliminary 
measurements  on  the  x-ray  reflectivity  from  the  surface  of  water.  These  measure¬ 
ments  were  then  continued  at  the  synchrotron  facility  in  July  of  1984. 

The  most  inportamt  accomplishment  of  the  preliminary  measurements  was 
the  development  of  sample  geometries  that  were  not  plagued  by  acoustically 
excited  capillary  waves.  These  results  were  published  in  the  Physical  BexJieij)  Lett. 
(Braslau  et  al. ,  Phys.  Eev.  Lett.  54,  114  (1985) .  Since  those  measurements  we  have 
been  continuing  to  explore  other  liquid  Surface  asd  have  made  extensive  measure¬ 
ments  on  CCl^  at  hASYLAB,  amd  methanol  at  Haurvaurd.  In  addition,  we  have  a 
further  measurement  on  the  follcwing  liquid  crystals:  9CB,  8S5  at  (HASYLAB); 

8CB  amd  QB7  (at  Harvard) ,  amd  one  microemulsion  system,  AOT-oil-water  as 
provided  by  John  Huang  of  Exxon.  The  CB7  sample  was  provided  by  Cyrus  Safinya 
of  Exxon,  and  he  paurticipated  in  the  measurement  at  Haurvaurd.  X-ray  reflectivity 
measurements  of  selected  solid  substrates  were  also  initiated  using  samples 
obtained  from  the  Harvard  Observatory.  These  samples  are  inportant  since  they 
are  protypical  of  surfaces  used  in  x-ray  mirrors.  We  have  seen  both  angular 
dependence  of  the  reflectivity  and  diffuse,  non-specrular  reflectivity.  We  are 
in  the  process  of  analyzing  this  data  in  order  to  determine  and  characterize 
the  roughness  of  solid  surfaces.  Cne  theoretical  model  has  successfully 
explained  most  of  the  observed  features;  however  there  cure  discrepancies  that 
indicate  the  model  is  not  complete.  X-ray  studies  of  freely  suspended  thin 
films,  and  of  the  Nematic-Smectic-A  phase  tremsition  have  continued,  however 
in  the  past  year  that  work  was  primarily  supported  by  other  grants  and  will 
not  be  discussed  further  in  this  report. 


1.13  Bt  Barriers  and  SNS  Proximity  Effect  Josephson  Devices.  O.  Liengme  ^ulcl 

M.  Tinkhan,  Contracts  N00014-84-K-0465  euid  N00014-83-K-0383;  Research 

unit  4  (former  #2) . 

Oonpaured  to  tunnel  junctions,  metallic  weak  links  present  negligible 
shunting  capacitemce ,  but  with  ordinary  metals,  they  usually  have  low 
resistamces,  which  make  ingpedance  matching  difficult  in  addition  to  causing 
severe  heating  problems  at  hig^  voltages  auid  frequencies.  It  has  been  pro¬ 
posed  that  a  aenrimetdl  barrier  might  conbine  some  of  the  advantages  of  metallic 
and  tunnel  barrier  devices.  In  experiments  now  in  progress,  we  find  that  Bi 
in  the  dirty  limit  is  well  described  by  the  standaurd  Ginzburg- Landau  theories 
used  for  ordinary  metal  SNS  junctions.  Ihis  holds  not  only  for  the  static 
properties,  such  as  the  temperature  dependence  of  the  critical  current,  but 
also  for  the  numerical  form  of  the  finite  voltage  I-V  curves,  as  computed 
with  the  TDGL  theory.  However,  these  experiments  also  show  that  the  low 
Fermi  velocity  chauracteristic  of  a  semimetal  with  low  carrier  density  leads 
to  a  short  coherence  length  unless  the  semimetal  is  very  cleam,  so  that  the 
Josefdison  coupling  is  usually  weak  for  liiiks  of  substauitial  length.  For 
realistic  finite  lengths,  mamy  questions  remadn  open.  For  example,  in  our 
simulations  auid  ejperiments  to  date  we  observe  a  cutoff  frequency  ~  2eI^I(/h. 
Hhereais  this  result  is  relatively  well  documented  in  the  case  of  short  devices , 
further  work,  both  eiperimental  and  analytical,  is  underway  to  complete  our 
understauiding  of  this  phenomenon  in  more  general  configurations  where  I^R 
falls  well  below  the  short  bridge  limiting  value. 
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II. I  Laser  Pulse  Shape  Storage  and  Reproduction.  N.W.  Carlson,  W.R.  Babbitt, 
Y.S.  Bai,  and  T.W.  Hossberg,  Contracts  N00014-84-K-0465  and 
DAAG-29-83-K0040;  Research  Unit  6  (former  #5). 

We  have  con^leted  a  series  of  experiments^  ^  which  demonstrates  that 
inhomogeneous ly  broadened  eUssorbers  can  "remenber"  the  temporal  shape  of  laser 
pulses  used  to  resonantly  irradiate  them.  We  have  also  shown  that  the 
absorbers  can  be  stimulated  to  emit  a  coherent  light  pulse  whose  shape  duplicates 
the  shape  of  an  earlier  excitation  pulse  (or  its  time  reverse ) .  Sairple-generated 
signals  have  been  observed  to  be  up  to  5  percent  as  intense  as  the  original  laser 
pulses.  If  one  considers  a  pulse  to  consist  of  a  digital  data  stream,  it  becomes 
obvious  that  the  cUssorber's  adiility  to  remember  to  pulse's  shape  may  be  important 
in  optical  data  storage  applications. 

In  our  first  experiment^  (which  is  basically  a  two-excitation-pulse  photon 
echo  experiment)  we  found  that  the  photon  echo  signal  constituted  a  time-reversed 
duplicate  of  the  first  laser  excitation  pulse  (pulse  1} .  We  performed  a  detedled 
study  of  the  conditions  necessary  to  ensure  pulseshape  reproduction.  We  found 
that  all  edisorbing  atoms  within  the  bemdwidth  of  pulse  1,  must  experience  pulse 
areas  less  than  or  approximately  equal  to  it/3>  Excellent  shape  reproduction 
was  found  when  the  total  duration  of  the  second  excitation  pulse  (pulse  2}  was 
less  them  the  temporal  modulation  time  of  pulse  1.  In  this  case,  the  spectral 
con^nents  of  pulse  2  are  flat  over  the  bandwidth  of  pulse  1.  Surprisingly,  we 
found  excellent  shape  reproduction  when  pulse  2  is  long  provided  that  it  is 


intense.  It  turns  out  that  shape  correlation  in  this  case  results  from  a  new 
type  of  power  broadening  in  which  cdjsorbers  at  various  frequencies  respond 
identically  in  terms  of  both  amplitude  emd  phase  to  a  strong  laser  field. 

Photon  echo  pulseshape  storage  is  mediated  by  am  optical  polarization 

which  unfortunately  means  that  it  camnot  persist  for  long  periods.  In  our 
2 

second  experiment,  however,  we  have  demonstrated  that  laser  pulseshapes  can 
be  stored  in  the  spectral  distribution  of  population  associated  with  a  single 
quantum  state.  Populations  may  be  extremely  staOsle  making  this  type  of  pulse- 
shape  storage  useful  in  optical  memory  application.  Storage  of  pulseshape 
information  in  populations  results  from  the  interference  between  two  laser 
pulses  separated  in  time.  The  pulseshape  information  is  recalled  later  by  a 
third  laser  pulse  which  stimulates  the  seunple  to  emit  a  duplicate  or  a  time- 
reversed  duplicate  of  one  of  the  first  two  laser  pulses.  This  processes  can 
be  viewed  as  a  generalization  of  the  stimulated  photon  echo.  Our  recalled 
signals  were  up  to  one  percent  as  intense  as  the  input  pulses. 

Finally,  we  demonstrated  that  laser  pulseshapes  can  be  stored  in  Zeemcm 
coherences.  The  storage  process  is  similar  to  that  described  in  the  previous 
paragraph,  but  the  first  two  pulses  are  made  to  excite  closely  spaced  adjacent 
tramsitions.  As  a  result,  pulseshape  information  is  stored  in  stable  coherences 
between  nearly  spaced  energy  levels. 

Additional  experiments  are  underway  to  elucidate  the  basic  mechemisms 
associated  with  the  above  processes,  and  to  demonstrate  that  they  will  be 
in^rtant  in  electro-optical  technology. 
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11.2  Real-Time  Optical  Waveform  Convolver/ Cross- Correlator.  Y.s.  Bai,  w.R. 
Babbitt,  N.W.  Carlson,  cuid  T.W.  Hossberg,  Contracts  N00014-84-K-0465 
and  DAAG29-83-K0040;  Research  Unit  6  (former  #5). 

We  have  performed  an  eiqseriment^  demonstrating  that  optical  coherent 
transient  effects  cem  be  employed  to  perform  convolutions  and  cross¬ 
correlations  between  the  tenoral  waveforms  of  optical  pulses.  Assume  that  a 
series  of  three  optical  pulses  reson^mtly  excite  inhomogeneous ly  broadened 
absorber.  If  one  of  the  pulses  is  made  teo^rally  short  while  the  other  two 
are  temporally  structured,  the  sanple  is  observed  to  emit  a  subsequent  pulse 
possessing  a  tezt^ral  envelope  proportional  to  the  square  of  the  convolution/ 
cross-correlation  of  the  electric  fields  of  the  ten^orally  structured  input 
pulses.  The  choice  between  convolution  jind  cross-correlation  is  determined 
by  which  of  the  three  input  pulses  is  made  teii?>orally  short.  Cptimized  output 
signals  were  about  one  percent  as  intense  as  the  input  signals.  In  addition 
to  our  experiments,  we  performed  numerical  simulations  to  determine  the  input- 
pulse  intensity  regimes  in  which  accurate  convolutions/cross-correlations  were 
obtained.  It  turns  out  that  accvirate  convolutions/cross-correlations  are 
obtained  only  when  atoms  in  the  absorber  respond  linearly  to  the  input  fields. 
This  basic  process  may  have  important  applications  in  optical  computing  systems. 
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Optical  Waveform  Oonvolver/Cross-correlator, "  Appl.  Phya.  Lett,  45,  714 
(1984). 

11.3  Optically  Inhibited  Collisional  Dephasing.  A.G.  Yodh,  J.E.  Golub,  N.W. 
Carlson,  cUid  T.W.  Mossberg,  Contracts  N00014-84-K-0465  and  NSF  PHY-82- 
0780;  Research  Unit  6  (former  #5). 

We  have  studied^  the  effect  of  a  strong  laser  field  on  gas-phase  colli¬ 
sional  relaxation.  We  find  that  a  strong  laser  field  acts  to  inhibit  the 
collisional  dephasing  associated  with  the  velocity-changing  aspect  of  collisions. 
The  laser  field  does  not  stop  collisions,  but  it  does  stop  some  of  the  phase 
randomization  normally  introduced  by  them.  Our  experiment  constitutes  the  first 
demonstration  that  gas-phase  elastic  collisional  relaxation  rates  can  be  laser- 
field  dependent.  Knowledge  that  the  collisional  dephasing  rates  are  laser-field 
dependent  may  be  important  in  the  modelling  of  gas  laser  dynamics. 
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11.4  Optical  Pulse  Compression  Using  Photon  Echoes.  Y.S.  Bai  and  T.W.  Mossberg, 
Contracts  N00014-84-K-0465  amd  DAAG29-83-K0040;  Research  Unit  6  (former  #5). 

We  have  performed^  a  theoretical  analysis  of  the  photon  echoes  which  may  be 
generated  with  frequency-chirped  optical  pulses.  We  find  that  two  linearly 
chirped  excitation  pulses  produce  echoes  whose  duration  is  roughly  equal  to  the 
inverse  chirp  bandwidth.  In  calculating  this  result,  we  assumed  that  both 
excitation  pulses  have  the  same  total  bandwidth,  but  that  the  chirp  rate  of  the 
second  pulse  is  twice  as  large  as  that  for  the  first  pulse.  The  compressed 
pulse  can  be  eiqiected  to  contain  on  the  order  of  ten  percent  as  much  energy  as 
was  contained  in  the  first  input  pulse.  The  echo  method  may  provide  a  very  use¬ 
ful  means  of  compressing  long  (i.e.,  microsecond)  optical  pulses. 
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11.5  Nonlinear  Four-Wave  Mixing  in  Vapors.  N.  Bloembergen  and  Y.H.  Zou, 

Contract  N00014-84-K-0465;  Research  Unit  7  (former  #6  and  part  of  #5). 

During  the  past  year  two  comprehensive  papers,  on  collision  induced 

coherence  emd  on  collision  induced  population  gratings  respectively,  have  been 
1  2 

published.  '  These  give  a  detailed  account  of  the  experimental  results  of 
four  wave  light  mixing  in  Na  vapor  with  various  buffer  gases.  The  frequency 
of  light  beams  is  detuned  by  about  1  cm  ^  from  the  D^  _  resonance  lines ,  so 
that  collision  induced  effects  are  clearly  dominant. 

During  the  past  year,  a  new  phenomenon,  the  collision  induced  Hanle 
effect,  has  been  observed  and  studied  in  some  detail.  This  effect  is  a 
variation  of  the  collision  induced  Zeeman  coherences  described  in  reference  1. 

A  simple  four-wave  mixing  configuration  is  used,  as  in  conventional  phase- 
conjugate  reflection.  All  light  beams  have  the  same  frequency  and  are  derived 
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from  a  single  dye  laser.  A  standing  wave  pattern  is  produced  by  two  light 
beams,  polarized  in  the  vertical  direction-  A  third  light  beam  makes  a  small 
angle  (about  0.3®)  with  the  first  two  and  is  polarized  in  the  horizontal 
direction.  The  beams  intersect  in  a  cell,  containing  Na  vapor  and  buffer  gas, 
placed  in  a  set  of  Helmholtz  coils.  The  intensity  of  the  phase-conjugate 
signal  polarized  horizontally  and  propagating  in  the  backboard  direction  of 
the  third  beam,  is  plotted  as  the  vertical  con^jonent  of  the  magnetic  field  is 
varied  through  zero.  Hanle-type  resonance  with  a  FWHM  width  of  about  30  milli- 
gauss  have  been  obtained.  The  peak  intensity  of  these  resonances  varies  as  the 
square  or  the  cube  of  the  buffer  gas  pressure.  These  first  results  on  the 
collision  induced  Hanle  effect  have  very  recently  been  published.^  A  systematic 
study  of  these  resonances  is  underway.  Their  parametric  dependence  on  the 
nature  and  the  pressure  of  the  buffer  gas,  the  detuning  from  the  resonance 

line,  the  power  level  and  polarization  in  the  three  incident  laser  beams,  the 
admixture  of  other  alkali-atoms,  will  be  investigated  further. 
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11.6  Two-Photon  Absorption  Spectroscopy  of  Rare  Earth  Ions.  c.D.  Cordero- 
Montalvo,  R.S.  Rana,  and  N.  Bloembergen,  Contract  N00014-84-K-0465 ; 
Research  Unit  7  (former  #6  and  part  of  #5). 

12  3+ 

During  the  past  year  the  level  of  the  4f  configuration  of  Pr 

in  LaCl-  has  been  found  for  the  first  time.^  This  forbidden  two-photon 

"*3  1  2 

transition  H^  S^,  with  Aj=4,  AL=  5  and  As=l,  has  also  been  studied 

for  the  same  ion  Pr^^  in  the  host  crystal  LaF  .  Furthermore  a  two-photon 

^  3+ 

absorption  study  of  both  cubic  non  non-cubic  sites  of  the  Gd  in  CaF^  has 
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been  convicted.  An  invited  review  paper  on  two-photon  spectroscopy  of  rare 
earth  ion  was  presented  at  the  International  Conference  on  Luminescence,  held 
in  Madison,  Wisconsin,  in  August  1984.  With  these  publications  the  activities 
in  this  particuleur  research  topic  have  been  terminated. 
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11.7  Picosecond  Laser  Interactions  with  Semiconductors.  H.  Kurz,  A.M. 

Malvezzi,  D.  von  der  Linde,  G.-Z,  Wang,  and  N.  Bloembergen,  Contracts 
N00014-84-K-0465  and  N00014-83-K-0030;  Research  Unit  7  (former  #6  emd  #5) . 

The  extensive  data  on  the  carrier  plasma  density  and  lattice  teiiperature 
during  and  following  picosecond  irradiation  of  silicon  sanples  has  been  reviewed 
and  analyzed  theoretically  in  eui  invited  paper^  at  the  meeting  of  the  Materials 
Research  Society,  held  in  Boston  in  November  1984.  A  theoretical  investigation 

of  the  variation  of  the  effective  mass  of  the  carriers  at  very  high  densities 

21  2 
(10  /cc)  and  ten^eratures  up  to  the  melting  point  has  also  been  carried  out. 

Experiments  on  the  heating  of  a  GaAs  crystal  under  picosecond  irradiation 
at  530  nm  wavelength  have  used  the  method  of  reflected  second  harmonic  genera¬ 
tion.  When  the  surface  melts,  the  production  of  second  harmonic  radiation  is 

3 

suddenly  reduced  .  This  effect  has  been  observed  and  reported.  It  was  con¬ 
cluded  that  the  dense  hot  carriers  exchange  energy  with  the  lattice  with  a 
characteristic  time  of  less  than  2  ps. 

Current  experiments  focus  on  the  evaporation  of  neutral  particles  from 

the  superheated  fluid.  It  has  been  established  that  the  silicon  surface 

2 

readily  reaches  a  temperature  over  5000®K  at  a  fluence  of  0.5  J/cm  at  530  nm 
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wavelength.  Thin  film  SOS  san^les  are  blcjwn  off.  This  blow-off  phenomenon 
has  hetmpered  the  study  of  the  rate  of  evaporation.  Experiments  are  now 
plemned  on  metal  films  of  aluminum,  and  possibly  silver  emd  gold.  The  lower 
melting  points  will  hopefully  permit  the  observation  of  a  range  of  super¬ 
heating  without  blow-off.  A  study  of  the  heating  and  evaporation  of  silicon 

4 

surfaces  coated  with  S1O2  layers  has  been  p\]blished. 

A  theoretical  study  was  also  made  to  obtain  shorter  pulses  (2-3  ps 
instead  of  20  ps)  from  a  Nd-Yag  mode-locked  laser  system  by  contoining  intra¬ 
cavity  self-phase-modulation  with  external  modulation.^ 
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and  Pulse  Con^ression  in  Mode- Locked  Lasers,"  Appt.  Phys,  A,  Springer- 
Heidelberg,  to  be  published,  1985, 


11.8  Picosecond  Laser- Induced  Emission  from  Metals  and  Semiconductors. 

H.  Kurz,  A.M.  Malvezzi,  and  N.  Bloembergen,  Contracts  N00014-84-K-0465 
auid  N00014-83-K-0030;  Research  Unit  7  (former  #6  and  #5). 

Careful  experimental  measurements  of  the  photo-electric  emission  induced 

by  X=  265  nm  picosecond  pulses  from  a  crystalline  silicon  surface  have  been 

con^ileted  over  a  wide  range  of  incident  fluences.  At  least  three  distinct 

2 

regimes  are  observed.  At  fluences  below  about  0.02  J/cm  ,  the  emitted  charge 
per  pulse  has  a  linear  suid  quadratic  component.  The  latter  is  due  to 


photo-emission  from  the  dense  hot  carrier  plasma,  llie  origin  of  this  nonlinear 

effect  was  confirmed  by  observation  with  temporally  separated  pulses.  The  synergetic 

effects,  when  the  two  pulses  overlap  in  time,  are  consistent  with  this  interpretation 

2 

In  the  second  regime,  with  fluences  in  the  range  0.02-0.15  J/cm  the  photo¬ 
emission  is  space  charge  limited.  This  is  in  agreement  with  a  theoretical 

2 

emalysis.  At  fluences  above  0.15  J/cm  ,  the  photo-emitted  electron  charge  in¬ 
creases  very  rapidly.  At  the  Scune  time,  a  strong  emission  of  positive  ions  is 
detected.  In  this  third  regime,  the  evaporation  during  the  picosecond  pulse  is 
sufficient  to  create  the  formation  of  a  plasma.  The  space  charge  effects  are 
neutralized  by  positive  ions.  The  velocity  of  these  positive  ions  has  also 
been  observed  by  a  time-of- flight  measurement.  There  is  an  indication  of  a 
fourth  regime,  in  which  ponderomotive  nonlinear  effects  in  the  created  plasma 
appear  to  become  dominant. 

These  investigations  have  been  presented^  at  the  Materials  Research 

Synposium  in  Boston,  November  1984.  A  detailed  account  has  also  been  accepted 

2 

for  pviblication  elsewhere. 

We  have  also  collaborated  with  the  group  of  Professor  Ippen  at  MIT  to 

study  the  photo-emission  by  75  femtosecond  pulses.  This  work  was  supported 

by  the  Joint  Services  Electronics  Program  at  MIT,  and  has  been  reported  in  a 

3 

recent  publication. 

With  these  reports  the  project  on  picosecond  photo- emission  has  been 
terminated. 
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11.9  Multi  photon  Vibrational  Excitation  of  Itolecules.  E-  Mazur,  J.  Wang, 
and  N.  Bloembergen,  Contract:  N00014-84-K-0465;  Research  Unit  8 
(former  #7). 

Time-resolved  Raman  spectroscopy  of  infrared  multiphoton  excited  mole¬ 
cules  has  proven  to  be  a  powerful  tool  in  the  study  of  intramoleculcu:  dyneunics. 
In  particuleu:,  information  on  the  intramoleculeu:  energy  distribution  can  be 
obtadned  as  a  function  of  time  and  infrared  laser  fluence. 

During  this  reporting  period  substauitial  modifications  and  improvements 
were  made  to  the  equipment  in  our  ladsoratory.  First,  the  high-power  picosecond 
CO^  laser  has  been  overhauled.  Ihis  was  necessary,  since  the  proposed  measure¬ 
ments  require  hundreds  of  thousemds  of  lasershots,  and  the  present  set-up  had 
been  operating  approximately  five  years  without  any  modification.  Second,  in 
order  to  svibstaintially  reduce  the  consumption  of  gas,  a  gas-recirculation 
system  is  being  added  to  the  CO2  laser  facility.  Finally,  the  old  computer 
system,  which  had  cui  inadequate  level  of  perf onocuice ,  has  been  replaced  by  a 
totally  new  computer-controlled  data-acquisition  system,  consisting  of  a 
Digital  computer  and  a  CAMAC  crate.  These  changes  constitute  major  improve¬ 
ments  which  will  make  better  and  more  efficient  measurements  possible. 

We  plan  to  continue  the  experimental  study  of  intramolecular  energy 
transfer  between  infrared  and  Raman  active  modes.  Measurements  will  be  per¬ 
formed  in  systems  where  several  Raman  active  modes  cem  be  observed,  in  order 
to  determine  if  the  excitation  of  these  modes  (after  infrared  multiphoton 
excitation,  but  before  amy  collisions  take  place)  corresponds  to  a  common 
temperature. 

1-3 

During  the  past  year  several  publications  have  either  appeared,  or 
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are  in  the  course  of  publication,  which  review  the  systematic  trends  of 
multiphoton  infrared  vibrational  excitation  as  a  function  of  molecular  size 
and  the  vibrational  density  of  states.  These  papers  are  based  on  the  large 
body  of  systematic  experimental  work  carried  out  under  this  project  over  the 
past  five  years. 
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1 1 1.1  Conditional  Density  Approach  to  Tracking.  R.W.  Brockett,  A.  Mclvor, 

and  K.  Wohn,  Oontracts  N00014-84-K-0465  euid  N00014-84“K-0504f  Research 

Unit  9  (former  #8) . 

Given  noisy  versions  of  successive  pictures  of  a  scene  which  conted.ns 
a  moving  cdsject,  it  is  natural  to  atten^t  to  confute  the  conditional  proba¬ 
bility  of  the  motion  peurameters.  These  parameters,  e.g.  the  tremslational 
cuid  rotational  velocities,  the  position  2tnd  orientation,  etc.,  may 
realistically  be  assumed  to  satisfy  stcdiastic  differential  equations.  The 
observations  vdiich  the  motion  generates  ccui  be  thought  of  as  the  gray  levels 
associated  with  the  individual  pixels  of  the  ccunera  which  tracks  the  motions. 
The  signals  which  underlie  these  observations  are  highly  correlated  and  the 
noise  tends  to  be  correlated  as  well.  Experiments  show  that  very  significeunt 
errors  result  if  these  correlations  are  ignored.  As  might  be  expected,  the 
differences  between  gray  levels  of  neighboring  pixels  gives  a  random  process 
with  much  less  spatial  correlation.  In  fact,  we  observe  that  it  is  often 
possible  to  get  satisfactory  results  by  modeling  these  differences  as  being 
independent  random  variables.  The  resulting  conditional  density  equation  has 
been  derived  and  propagated  in  some  simple  situations.  This  pr<^abilistic 
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model  seems  to  work  well  cuid  extensions  to  more  con^lex  problems  axe  being 
pursued. 

1 1 1 .2  Probabilistic  Methods  in  Single  Frame  Picture  Processing.  R.w.  Brockett, 
A.  McXvor,  and  W.  Jian,  Contracts  N00014-84-K-0465  euid  NQ0014-84-K-0504; 
Resecurch  Unit  9  (former  #8) . 

The  work  of  Alan  Mclvor^  is  devoted  to  a  two-dimensional  innovations 
approach  to  picture  processing.  It  is  rather  intuitive  that  if  one  has  a 
jump  process  which  switches  between  two  Gauss  Meurkcv  processes  then  a  whitening 
filter  followed  by  a  demodulator  would  be  a  suitable  "perfect"  detector  of  the 
jump  process.  We  have  shown  that  veurious  algorithms  in  use  in  picture  pro¬ 
cessing  ceui  be  interpreted  as  approximations  to  this. 

We  have  also  investigated  probabilistic  approaches  to  edge  tracking. 

This  involves  constructing  a  probabilistic  model  for  edges.  In  this  area  we 
have  investigated  two  possibilities,  one  beised  on  modeling  the  curvature  as 
a  Gauss  Markov  process  and  the  other  based  on  modeling  it  as  a  finite  state 
process.  Results  of  con^uter  simulations  indicate  that  while  there  is  some¬ 
thing  to  be  gained  from  each  approach,  sharp  corners  cause  problems.  We  have 
also  investigated  a  two  point  boundeury  value  problem  type  interpolation  scheme 
to  continue  edges  through  faint  regions.  This  work  will  be  reported  on  in  the 
forthcoming  thesis.^ 

Reference 

1.  Alan  Mclvor,  Probabilistic  Models  in  Computer  Visions,  Ph.D.  Thesis, 
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NI.3  Mul ti mode- Non 1  inear  Control.  R.w.  Brockett,  Contracts  N00014-84-K-0465, 

DAAG29-79-C-0147,  AFSOR-81-7401,  ECS-81-21428 j  Research  Unit  9  (former  #8) 

Recent  work  on  nonlinear  control  has  emphasized  problems  for  which  the 
nonlinear  effects  are  amalytic  functions.  On  the  other  hwd,  static  friction, 
one  of  the  most  frequently  encountered  nonlinear  effects ,  is  definitely  not 
analytic.  In  a  recent  paper^  we  undertook  a  study  of  what  we  call  "multimode" 
systems— a  class  which  includes  circvdts  with  diodes  and  gear  trains  with 
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deadzone,  backlash  emd  friction.  Ihe  main  results  in  this  paper  give  a 
measure  of  how  much  performance  degradation  such  effects  cause.  In  addition, 
we  present  a  convenient  formalism  for  describing  such  systems.  The  most 
tractable  class  of  useful  models  seems  to  be  that  of  piecewise  lineau:  systems 
with  the  added  feature  that  in  some  "boundary  regions”  the  order  of  the 
system  drops.  Ihis  reflects,  for  example,  a  diode  shorting  out  a  capacitor 
when  its  voltage  goes  positive  or  gears  leaving  a  backlash  region.  Ihis 
formalism  suggests  a  convenient  simulation  description  as  well. 


Reference 
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1984. 


III. 4  Discrete  Event  Dynamic  Systems  Study.  Y.C.  Ho,  R.  Suri,  X.  Cao,  G. 
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We  continue  the  study  described  last  year  at  the  triannual  JSEP  review 
on  discrete  event  systems.  Specifically: 

In  the  area  of  perturbation  analysis  (P/A)  of  13ED6,  we  resolved  certain 

open  theoretical  question  concerning  the  interchamge  of  differentiation  emd 

expectation  operators  for  Vi/G/1  queueing  systems.  We  also  showed  that  P/A 
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could  be  used  for  very  fast  optimization  of  m/G/1  queues.  We  have  also 
presented  an  analytical  proof  that  P/A  is  much  more  efficient  than  alter- 

3 

native  techniques  for  sensitivity  walysis  of  a  general  class  of  systems. 

In  the  curea  of  manufacturing  systems,  we  have  continued  our  efforts  in 
the  study  of  flexible  numufacturing  as  well  as  providing  system-theoretic 

9 

approaches  to  manufacturing  systems  in  general. 

12 

In  the  communication  network  area,  we  conpleted  a  study  of  decen¬ 
tralized  flow  control  based  on  the  concept  of  a  market  mechemism  for  the 
allocation  of  scarce  resources.  The  basic  idea  is  to  make  each  user  decl£u:e 
what  he/she  needs  2uid  to  allocate  optimally  based  on  these  decleurations .  The 
important  things  here  is  that  the  payment  scheme  for  the  resources  is  designed 


to  induce  each  user  to  report  truthfully  in  his/her  own  interest.  Further¬ 
more,  no  central  authority  is  needed  to  make  such  allocations.  Local  flow 
control  decision  based  on  delayed  reporting  of  declaration  from  other  sources 
cam  be  used  and  proved  to  yield  stable  allocations  over  time. 
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Resecurch  in  the  area  of  electromagnetic  radiation  is  directed  toward  the 
solution  of  practical  problems  through  the  con^lete  understemding  of  the  under¬ 
lying  physical  phenomena.  This  involves  the  coordinated  application  of  modern 
analytical,  numerical,  and  experimental  techniques  and  the  use  of  high-speed 
computers  emd  precision  instrumentation.  Application  is  also  made  of  modeling 
techniques  and  the  principle  of  similitude.  Most  practically  significant  prob¬ 
lems  in  the  area  are  sufficiently  complicated  that  extensive  computation  emd 
measurement  are  often  required  to  justify  approximations  that  are  usually 
necessary.  Where  possible,  general  formulas  are  obtained  and  verified  experi¬ 
mentally  so  that  the  phenomenon  under  study  can  be  understood  physically  in 
analytical  form  and  not  just  as  a  set  of  numbers. 

The  researches  eure  concerned  primarily  with  the  properties  of  antennais  and 
arrays  amd  of  the  electromagnetic  fields  they  generate  in  various  practically 
importamt  environments  that  lead  to  difficult  problems  with  complicated  boundary 
conditions .  Examples  include  dipoles ,  traveling-wave  antennas  and  arrays , 
crossed  dipoles,  amd  loops  near  the  boundary  between  two  media  such  as  air  amd 
eaurth  or  sea  water,  and  rock  and  sea  water;  the  scattering  of  electromagnetic 
waves  from  buried  or  svibmerged  objects ;  the  properties  of  lateral  electromag¬ 
netic  waves  and  pulses;  the  field  inside  and  currents  and  chaurges  induced  on 
obstacles  within  an  emp  simulator;  arrays  of  antennas  along  curves  lines;  and 
solitary  electromagnetic  pulses  with  slow  decay. 
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IV. 1  Theoretical  Study  of  Scattering  from  Buried  or  Submerged  Objects.  R.w.P. 
King,  T.T.  Wu,  and  W.-Y-  Pein,  Contract  N00014-K-0465;  Research  Unit  11 
(former  #10) . 


The  detection  amd  localization  of  buried  or  submerged  objects  or  inhomo¬ 
geneities  in  a  dissipative  half-space  from  measurements  of  a  scattered  electro¬ 
magnetic  field  are  of  current  interest.  Earlier  investigations  have  been 
concerned  primarily  with  plane  waves  incident  on  the  surface  of  the  earth  from 
the  air  above  at  an  arbitrary  amgle.  The  present  study  deals  with  the  incident 
fields  generated  by  horizontal  and  vertical  antennas  located  on  the  surface  of 
the  earth.  The  first  step  in  the  analysis  is  the  determination  of  this  incident 
field.  The  complete  electromagnetic  field  paurallel  to  the  surface  of  salt  water 
has  been  determined  '  «hen  the  source  is  a  horizontal  electric  dipole  either 
with  open  ends  or  grounded  at  its  ends.  The  dipole  is  in  the  zdr  at  a  small 
height  d  above  the  water.  Actually  such  an  amtenna  is  equivalent  to  am  eccen¬ 
trically  insulated  amtenna  lying  on  the  surface  of  the  half-space.  The  effective 
thickness  of  the  insulation  is  d.  This  type  of  antenna  has  the  properties  of  a 
transmission  line  with  a  coitftlex  wave  number  and  characteristic  impedamce.  By 
determining  the  distribution  of  current  and  the  effective  length  of  the  amtenna, 
the  entire  field  along  the  surface  of  the  earth  can  be  readily  obtained  with  the 
help  of  newly  derived  comprehensive  formulas.  The  components  of  primary  interest 
aire  those  tamgent  to  the  boundary  surface,  viz.,  E^p  or  and  E^^  or  B^p. 

The  field  of  the  open-ended  dipole  is  found  to  be  simpler  than  that  of  the  ter¬ 
minated  dipole  since  the  currents  in  the  vertical  terminations  of  the  latter 
generate  a  significauit  field  consisting  of  the  three  con^onents  Ep,  amd  B^. 
The  combined  field  of  the  horizontal  and  vertical  members  of  the  coi[4>lete  antenna 
can  be  quite  complicated  if  the  height  d,  and  with  it  the  length  d  of  the 
vertical  elements ,  is  sufficient  to  make  the  contributions  from  them  coit^arable 
to  that  from  the  horizontal  element.  At  close  range  Ep  and  E^  are  comparable 
in  magnitude.  At  distances  for  which  k2P>l,  Ep  is  significantly  greater 
since  it  decreaises  only  ais  1/p  whereas  E^  decreases  as  1/p^,  The  cylindrical 
con^Jonents  Ep  and  E^  aure  individually  simpler  than  the  related  Cartesian 
components  E^  and  E^. 

Once  the  incident  field  is  known,  the  scattered  field  can  be  investigated. 
The  first  scattering  object  studied  is  the  horizontal  insulated  conductor.^  This 
provides  a  scattered  field  that  is  readily  calculated  and  combined  with  the 


incident  field  to  obtain  a  total  field  that  exhibits  interference  phenomena. 

Ihe  length  of  the  insulated  metal  rod  was  chosen  to  be  near  resonance  at  the 
operating  frequency  for  the  r2uige  of  submerged  depth  d  considered.  The  inci¬ 
dent  lateral-wave  field  is  maintained  by  a  horizontal  dipole  a  height  d^  in 
edr  above  the  surface  of  the  earth,  as  described  above.  The  total  field,  i.e., 
the  incident  plus  scattered  fields,  is  observed  in  the  air  also  just  above  the 
surface  of  the  eeurth.  The  current  Induced  in  the  insulated  wire  by  the  lateral- 
wave  field  incident  from  the  transmitter  is  evaluated  first  and  then  the  re¬ 
radiated  field  just  above  the  surface  is  calculated.  Two  orientations  are 
considered:  The  field  at  the  center  of  the  buried  conductor  is  (1)  parallel  to 
Eip(p,0,z),  emd  (2)  parallel  to  (P,tt/2,z)  .  The  incident,  scattered  and 

total  fields  in  an  area  above  the  conductor  have  been  evaluated  for  the  condi¬ 
tions  of  a  laboratory  model.  It  is  clearly  demonstrated  that  a  significant 
change  in  the  electric  field  over  the  volume  occupied  by  the  scattering  object 
is  observed  and  that  this  is  sufficiently  localized  to  permit  an  accurate 
bounding  of  the  area  above  the  object.  The  scattering  from  objects  like  a  metal 
disk  or  any  region  with  a  wave  number  that  is  significiintly  different  from  that 
of  the  ambient  medium  will  be  investigated  quantitatively  since  the  method  has 
definite  promise  for  discovering  and  locating  underground  oil  fields.  Extensive 
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measurements  in  a  salt-water-filled  tank  containing  a  submerged  metal  disk  have 
demonstrated  the  feasibility  of  the  method. 

Only  a  few  scattering  problems  involving  very  simply  shaped  targets  can 
be  solved  analytically.  In  a  separate  study  by  W.-Y.  Pan,  analytical  expressions 
have  been  obtained  for  the  scattered  field  from  a  buried,  perfectly  conducting 
rectangular  plate  placed  parallel  to  the  interface  between  air  and  earth. ^  First 
the  current  is  expressed  in  terms  of  several  Tchebichef  polynomials  of  the  second 
kind.  The  scattered  field  is  the  s\am  of  the  fields  that  are  generated  separately 
by  equivalent  horizontal  dipoles  and  multipoles.  The  moments  of  the  dipoles  and 
multipoles  are  obtained  from  the  coefficients  of  the  current  expansion.  These 
coefficients  have  been  calculated  over  a  frequency  range  from  13  to  160  kHz  and 
a  range  of  plate  sizes  from  2  to  20,000  square  meters  in  area. 
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Because  the  recently  derived  sets  of  formulate  for  the  conplete  electromag- 

1-3  4 

netic  fields  of  horizontal  auid  vertical  electric  dipoles  on  or  near  the  bound¬ 
ary  between  two  half-spaces  are  both  simple  and  accurate,  they  have  provided  new 
insight  into  lateral-wave  propagation  in  general  and  more  flexible  tools  for  a 
variety  of  specific  applications.  Examples  of  their  usefulness  are  presented  in 
several  new  papers.  One  paper^  describes 'a  number  of  practical  antennas  poten¬ 
tially  useful  in  communicating  with  submarines.  These  include  bare  vertical  mono¬ 
poles  cUid  arrays  of  such  monopoles  in  air  on  the  earth's  surface;  horizontal 
traveling-wave  emtennas  (Beverage  antennas)  or  arrays  of  such  antennas  in  air 
close  to  the  surface  of  the  earth  or  the  sea;  and  insulated  traveling-wave  ele¬ 
ments  in  the  sea  used  either  singly  or  in  arrays  designed  to  provide  a  radially 
progressive  field  with  a  maximum  vertically  upward.  The  properties  of  the 
antennas  are  evaluated  in  the  frequency  range  10 < f <  30  kHz  which  is  optimum 
for  depths  up  to  50  m. 

A  second  paper^  demonstrates  even  more  clearly  that  the  new  formulas  for 
the  surface-wave  fields  of  antennas  near  a  boundary  surface  have  provided  an 


expanded  horizon  for  understanding  and  using  surface  waves.  By  examining  the 
formulas,  surface  waves  are  seen  to  have  unusual  properties  that  make  them 
valuable  tools  in  geophysical  prospecting.  The  formulas  show  that  over  all 
distances  and  at  all  frequencies  a  signal  reaching  a  receiving  dipole  in  the 
salt  water  near  the  sea  floor  has,  in  effect,  proceeded  downward  from  the  trems- 
mitting  dipole  in  the  sea  water  into  the  rock,  then  radially  outward  in  the 
rock,  and  finally  upward  and  back  into  the  sea  water  to  the  receiving  dipole. 

The  decrease  in  amplitude  and  chamge  in  pheise  with  radial  distance  depend  prima¬ 
rily  on  the  electrical  properties  of  the  rook.  With  a  suitable  choice  of  fre¬ 
quencies  and  radial  distances  p  between  transmitting  emd  receiving  dipoles, 
the  received  signal  cam  have  the  simple  form  exp(-ap)e}q>(iBp)/P,  where  a  is 
the  attenuation  constant  and  B  is  the  real  wave  number  of  the  rock.  Ohese 
important  facts  are  the  basis  for  the  development  of  methods  for  the  determina¬ 
tion  of  the  conductivity  and  permittivity  of  the  rock  forming  the  sea  floor  by 
means  of  measurements  made  in  the  adjacent  sea  water.  Since  lateral  waves  do  not 
penetrate  deeply  into  the  ocean  floor,  they  au:e  useful  primarily  in  determining 
the  average  conductivity  and  permittivity  of  the  materail  fairly  close  to  the 
boundary.  For  this,  a  two- layer  theory  may  be  adequate.  In  this  same  paper 
the  new  formulas  are  used  to  assist  in  the  interpretation  of  available  measure¬ 
ments  of  the  conductivity  of  the  sea  floor.  It  appears  that  the  measured  data 
would  be  well  satisfied  by  a  half-space  model  of  the  lithosphere  that  is  homo¬ 
geneous  but  one-dimensionally  anisotropic  in  the  conductivity.  Such  a  model 
involves  no  reflections  from  deeper,  widely  spaced  horizontal  layers.  An  ana¬ 
lytical  study  by  W.-Y.  Pan  of  the  one-dimensionally  anisotropic  half-space  is 
now  in  progress.  Preliminary  results  show  that  with  a  a  “  compo- 

nents  and  of  electric  type  depend  primarily  on  a^,  while  the 

coiq}onents  magnetic  type  depend  primarily  on  a^  =  a^. 

Ihis  behavior  agrees  with  the  available  measured  data.  A  paper  on  this  subject 
is  being  prepared. 

A  summary  paper  which  discusses  these  various  applications  of  lateral 
electromagnetic  waves  has  been  presented  at  a  recent  scientific  meeting^  emd 
submitted  for  a  special  issue  of  the  IEEE  Transactions  on  Antennas  and  Propaga- 

Q 

tion. 

Currently  the  Wu-King  formulas  are  being  used  to  determine  certain  prop¬ 
erties  of  lateral  waves ,  such  as  the  maximum  depth  of  penetration  into  Region  2 


of  the  lateral-wave  part  of  the  field,  auid  the  fraction  of  power  associated 
with  the  lateral  wave  in  a  specified  range.  To  accoit^jlish  this,  it  is 
necesseury  to  obtain  e:q>licit  formulas  for  the  complete  field  in  Segion  2 
and,  from  them,  the  locus  of  the  Poynting  vector.  These  have  been  derived, 
and  a  paper  applying  this  information  to  the  use  of  horizontal  2uitennas  on 
the  sea  floor  for  geophysical  exploration  of  the  lithosphere  is  being  written. 
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The  propagation  of  electromagnetic  waves  along  boiuidaries  between  quite 
different  media  is  predominamtly  by  lateral  waves.  This  is  well  understood 
along  smooth  boundcuries ,  but  the  scattering  of  lateral  waves  at  boundaries 
and  discontinuities  is  not  well  understood.  A  theoretical  study  has  been 
made^  of  the  effect  on  the  field  generated  by  a  horizontal  dipole  due  to 
vertical  discontinuities  on  both  sides  of  the  boundcury  between  earth  (water) 
and  ed.r.  The  two  types  of  reflecting  boundaries  considered  are ;  a)  vertical 
discontinuities  in  the  denser  medium  (earth,  water)  in  the  form  of  changes  in 
the  wave  number  from  to  k^;  and  b)  vertical  projections  into  the  less 

dense  medium  (air)  in  the  form  of  thin  walls  of  the  denser  medium.  The  analysis 
shows  that  reflections  from  the  discontinuities  in  the  denser  half-space  are 
extremely  small  even  when  k^  changes  from  a  value  cheuracteristic  of  salt 
water  to  a  value  (k^ •*•“’)  characteristic  of  a  perfect  conductor  (metal). 
Reflections  from  metal  walls  in  air  when  *3^®  much  larger  but  still 

quite  small  vhen  the  height  of  the  wall  is  a  tenth  of  a  wavelength  or  less. 

They  are  very  significant  when  the  wall  is  a  quarter-wave  length  hi^.  These 
results  cure  readily  understood  if  it  is  recalled  that  the  vertical  unit  dipole 
in  air  is  very  superior  to  the  vertical  unit  dipole  in  salt  water  as  a  generator 
of  lateral  waves.  Steuiding-wave  patterns  due  to  incident  and  reflected  lateral 
waves  cure  determined  emd  shown  graphically. 

The  reflections  of  lateral  waves  at  cylindrical  discontinuities  above  emd 
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below  the  salt-water  surface  have  been  measured  in  a  semicirculeu:  tank  at  two 
frequencies,  f » 644  MHz  euid  f^l.S  GHz.  Single  combinations  of  horizontal 
and  vertical  metal  pleuies  were  chosen  as  the  discontinuities.  The  measurements 
confirm  the  theoretical  predictions  discussed  above.  The  small  reflections 
from  the  ceiling  were  also  studied  e:q)erimentally  and  theoretically  in  their 
effect  on  both  traveling  and  standing  lateral  waves.  The  meeisured  field 

patterns  involve,  in  general,  three  waves:  1)  the  direct  lateral  waves  with 
both  vertical  and  radial  consonants  traveling  radially  outward  from  the  source; 
2)  reflections  from  the  cylindrical  boundary  in  the  form  of  lateral  waves 
traveling  radially  inwaurd  toward  the  source  and  producing  a  reguleu:  standing- 
wave  pattern;  and  3)  transverse  plane  waves  coming  from  the  ceiling  at  normal 


or  near  normal  incidence  £md  polarized  with  the  electric  field  in  the  radial 
direction.  These  produce  interference  patterns  with  traveling  and  steuiding 
lateral  waves  that  are  quite  different  from  ordinary  standing  waves  produced 
by  simileu:  waves  traveling  in  opposite  directions.  It  appears  that  the 
con^lications  introduced  by  the  finite  size  6uid  indoor  location  of  the  tank 
are  now  well  understood. 

In  a  separate  but  related  study,  the  lateral -wave  field  is  incident 
upon  6md  scattered  by  a  thin  metal  cylinder.  The  usefulness  of  the  new  euia- 
lytical  formula  for  the  z-con^onent  of  the  electric  field  and  the  generaliza¬ 
tion  of  the  theojry  to  fatter  cylinders  are  discussed  in  a  recently  presented 
paper. ^  The  measurements  of  both  incident  and  scattered  fields  for  thin  and 
fat  cylinders  near  the  interface  between  two  regions  of  a  model  lithosphere, 
one  or  both  of  which  may  be  described  by  complex  wave  numbers,  are  presented. 

References 

1.  R.W.P.  King,  "On  the  Reflections  of  Lateral  Electromagnetic  Waves  from 
Perpendicular  Boundeuries , "  tf.  of  Appt.  Phys.  5S ,  3916-3926  (June  1984). 

2.  M.F.  Brown,  B.W.P.  King,  auid  T.T.  Wu,  "Experiments  on  the  Reflection  of 
Lateral  Electromagnetic  Waves,"  J.  of  Appl.  Phya.  55,  3927-3933  (June  1984). 

3.  M.F.  Brown,  "Scattering  of  a  Lateral  Electromagnetic  Wave  by  a  Metal 
cylinder,"  presented  at  the  National  Radio  Science  Meeting,  Boston,  MA, 

June  1984. 


IV. 4  Lateral "Wave  Propagation  and  Modeling  of  the  Lithosphere.  M.F.  Brown, 
R.W.P.  King,  T.T.  Wu,  and  J.T.  deBettencourt ,  Contracts  N00014-79-C-0419 
and  N00014-84-K-0465;  Reseeurch  Unit  11  (former  #10). 

The  first  of  a  series  of  experiments  which  will  culminate  in  an  experi¬ 
mental  and  theoretical  treatment  of  lateral-wave  pulses  has  been  undertaken. 

The  prelimin5u:y  problem  of  determining  the  fundamental  tremsient-response 
characteristics  of  dipole  antenncis  in  an  indoor  environment  has  been  studied. 

In  addressing  this  problem,  two  vertical  monopoles  have  been  used  above  a 
metal  ground  plane.  A  receiving  monopole  has  been  designed  to  function  in  a 
mode  whidi  eliminates  noise  by  utilizing  a  twin-measurement  procediire  developed 
in  conjunction  with  emp  simulator  measurements.^  In  such  a  procedure,  the 


received  signal  is  measured  twice,  once  with  the  monopole  pointing  "down" 

(at  the  ground  plcute)  and  once  "up".  The  signals  are  processed  by  a  waveform 
computer  programmed  to  subtract  the  two  signals  at  its  input  so  that  amy 
induced  currents  in  the  laboratory  which  axe  not  part  of  the  transmitted  pulse 
are  eliminated. 

Ihe  pulse  waveforms  of  the  received  signal  for  v£u:ious  separations  between 
the  dipoles  have  been  plotted  together  with  the  aissociated  Fourier  transforms. 
Plots  of  the  electric  field  versus  time  have  proved  to  be  quite  conplicated, 
i.e.,  with  multiple  maxima  and  minima,  for  the  first  choice  of  tramsmitting- 
cuitenna  length  (13.5  cm).  The  Fovurier  transform  of  this  field  shows  a  maximum 
at  ~  500  MHz,  while  the  shortest  wavelength  conponent  is  ~  .8  GHz-  The 
origins  of  the  cosplicated  pulse  sequence  have  been  investigated  by  studying 
i)  the  resoneuice  qualities  of  the  dipole  tremsmitting  antenna,  ii)  the  effect 
of  a  vertical  ground  plcuie  (image  pleuie)  close  to  the  dipole,  andi  iii)  the 
length  of  the  ouitenna  as  it  relates  to  the  width  of  the  pulse  (pulse  at  half¬ 
amplitude  is  'xlns).  The  latter  two  have  been  found  to  be  the  most  important 
in  influencing  the  pulse  shape.  The  complications  in  this  pulse  shape  have 
been  considerably  reduced  by  significantly  lengthening  the  transmitting  emtenna. 
A  distinct  peak  pulse  may  now  be  resolved;  smaller  undulations  corresponding  to 
reflections  having  long  delay  times  (relative  to  that  associated  with  the  inter- 
cmtenna  separation)  are  still  present. 

Current  work  is  being  directed  toward  further  refining  and  sioplifying  the 
received  waveform  by  moving  the  source  dipole  closer  to  the  horizontal  center 
of  the  underlying  metal  surface  assembled  in  the  laboratory.  In  this  way  it  is 
expected  to  be  less  influenced  by  reflections  from  the  vertical  ground  plane 
neeu:  which  it  lies.  Such  reflections  axe  quite  complicated  since,  although 
image  theory  applies ,  an  image  is  associated  with  each  frequency  component  of 
the  radiated  signal. 

As  a  fxirther  check  on  the  apparatus,  the  vertical-dipole  transmitting 
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element  has  been  connected  to  the  vertical  groiuid  plane  of  earlier  studies  ' 
near  which  near  which  it  lies,  so  that  the  euitenna  forms  essentially  a  short 
circuit  as  opposed  to  an  open  circuit.  In  this  way  the  pulse  is  reflected 
back  with  opposite  poleurity  to  that  observed  for  the  ungrounded  ceise  6uid 
produces  a  noticeable  difference  in  the  leading  edge  of  the  first  meocimum  in 
the  pulse  waveform. 


The  second  set  of  measurements  in  the  aforementioned  sequence  is  about 
to  commence,  naunely,  the  study  of  vertical  monopoles  above  salt  water.  This 
represents  a  transition  to  a  half-space  problem.  The  waveforms  are  expected 
to  be  quite  similar  to  those  obtained  over  the  metal  underburden  (ground 
plane)  discussed  above.  The  ground-plane  studies  have  thus  been  made  to 
serve  two  purposes:  i)  to  test  the  apparatus,  and  ii)  to  serve  as  a  basis 
for  comparison  with  the  salt-water  results. 

In  a  separate  study  an  investigation  has  been  made  of  the  effects  of 
atmospheric  noise  on  the  lateral-wave  propagation  between  submerged  horizontal 
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dipoles  in  the  sea.  Calculations  have  been  made  for  sea  water  with  ~ 

0j^=  4  S/m.  The  vertical  noise  is  for  the  worst-time  conditions  of  summer-time 
for  Eastern  U.S.A.  at  40'’N  latitude.  Curves  are  shown  for  the  vertical  emd 
horizontal  noise  components  of  the  electric  field  at  the  surface,  and  the  field 
strengths  of  the  radial  component  of  the  lateral-wave  field  generated 

by  an  infinitesimal  horizontal  electric  dipole  with  unit  moment  in  the  sea  at 
radial  distcuices  from  p=l  to  100  km.  Similar  curves  cure  also  presented  for  a 
practical  anteaina  system  with  moment  lAS.  =»  100  x  100  *  10^  A*m.  The  resulting 
signal-to-noise  ratios  have  been  plotted  as  functions  of  the  radial  range  for 
f-1,  10,  100  and  1000  kHz.  These  results  infer  useful  communication  to  sub¬ 
marines  using  lateral  waves  limited  by  attenuated  atmospheric  noise  for  shallow 
depths.  For  greater  depths,  lower  frequencies  should  provide  more  optimxam 
signal-to-noise  ratios. 
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Current  interest  in  the  possibility  of  investigating  the  electrical 
properties  of  the  earth  beneath  the  oceeui  by  electromagnetic  mezms  naturally 
leads  to  the  problem  of  wave  propagation  in  layered  media.  A  realistic  model 
of  the  earth's  crust  typically  includes  a  layer  of  sediment  followed  by  mzmy 
layers  of  rock  of  various  types.  Lateral-wave  propagation  for  such  czises  czui 
be  much  more  complicated  since  a  different  wave  cam  travel  in  each  of  the  rock 
regions  and,  at  the  same  time,  each  of  the  waves  is  influenced  by  the  neigh¬ 
boring  regions.  Geophysicists  are  interested  in  the  electrical  properties  of 
these  rock  layers  as  well  as  their  thicknesses  and  orientations. 

Hie  simplest  multi-layered  model  that  takes  into  account  the  presence  of 
a  hi^ly  conducting  mantle  region  beneath  a  layer  of  rock  is  that  of  a  hori¬ 
zontal  dipole  zuitenna  located  near  a  salt  water-air  interface  with  a  perfectly 
conducting  top  plate  above  the  air  region.  Hie  salt  water  represents  the 
ocean,  the  air  region  is  the  analog  of  the  rock  layer,  and  the  top  plate 

corresponds  to  the  highly  conducting  mantle.  A  theoretical  investigation  of 

1-3 

this  problem  has  been  czurried  out.  Hie  theoretical  expression  for  the 

radial  electric  field  is  in  the  form  of  modes  and  may  be  obtained  either 

P 

directly  from  the  integrals  or  from  approximate  perturbative  methods.  It  is 
no  coincidence  that  the  expression  for  is  very  similar  to  that  for  a 

waveguide  since  the  aU.r  region  is  very  much  like  a  slab  waveguide  with  a 
lossy  wall  at  the  air-water  interface. 

An  experimental  study  of  lateral-wave  propagation  in  the  presence  of  a 
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top  plate  has  also  been  completed.  '  '  As  predicted  by  theory,  the  measured 
fields  eidiibit  a  waveguide- type  behavior,  i.e.,  the  pattern  is  understandable 
in  terms  of  modes.  For  the  top  plate  at  a  height  h  above  the  water-filled 
tzmk  of  less  than  a  half-  .velength,  an  e3q>onential  decay  is  observed  with 
small  superimposed  oscillations  due  to  side-wall  reflections.  For  heights 
between  a  half  and  a  full  wavelength,  another  mode  is  observed  which  produces 
large  oscillations  in  the  pattern.  As  the  height  is  increased  further,  more 
and  more  modes  come  into  play,  zis  expected. 


A  more  realistic  three-layer  model  consists  of  a  semi-infinite  water 
layer  (Region  1) ,  a  sediment  layer  composed  of  a  mixture  of  samd  aind  water 
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(Begion  2} ,  euid  a  semi-infinite  rock  layer  (Begion  3) .  A  theoretical  study 
of  this  configuration  hats  been  made.  '  '  An  infinitesimal  horizontal  time- 
harmonic  electric  dipole  is  in  Begion  1  a  distaince  d  above  the  interface 
between  Begions  1  and  2  (z^O).  Analytical  expressions  have  been  obtained 

for  the  electromagnetic  fields  in  Begion  1  near  the  interface,  subject  to  the 
following  restrictions:  (1)  Both  d  amd  z  axe  small  conpaired  to  a  wavelength 
in  the  water,  where  d  is  the  height  of  the  tramsmitting  dipole  above  the 
interface  (z^O)  between  Begions  1  amd  2  and  z  is  the  height  of  the  receiver 

(2)  the  radial  distance  p  from  the  dipole  ic  large  compared  to  z  and  d;  auid 

(3)  the  two  parameters  a^k^/k^  and  £=k^/k2  are  small  in  magnitude,  where 
k^  is  the  coitplex  wave  number  of  the  j-th  medium.  Ihese  three  criteria  aure 
satisfied  nicely  in  practical  causes.  The  three- layer  solution  looks  very  much 
like  the  Wu-King  two- layer  solution.  The  differences  aure:  (1)  there  is  a 
sediment  term  which  represents  the  contribution  from  the  wave  traveling  in  the 
sediment  layer;  at  distauices  and  frequencies  of  interest,  it  is  negligible; 

(2)  there  is  a  constauit  A  which  cam  chamge  the  magnitudes  of  the  fields;  amd 

(3)  the  Fresnel  integral  now  depends  on  a  complicated  constant  V.  It  is  this 
greater  conpli cation  that  makes  the  solution  interesting.  Hopefully  a  meams 
cam  be  determined  to  exploit  these  differences  amd  ther^y  determine  not  only 
the  paurauneters  of  the  rock  region  but  also  the  thickness  of  the  sediment  layer. 

current  work  on  this  project  is  being  devoted  to  modifying  the  experimental 
setup  used  in  the  top-plate  measurements  so  that  the  more  importamt  case  of  adr 
over  a  sheet  of  dielectric  over  water  can  be  studied.  All  of  the  necessau:y 
equipment  is  avadlable  including  the  model  tank,  generator,  receiver,  amd  auto¬ 
mated  movable  receiving  pr<^e.  Construction  of  an  intermediate  layer  using 
sheets  of  wallboaurd  of  different  thicknesses  is  required.  To  begin  with  both 
the  tramsmitting  amd  receiving  amtennas  will  be  located  in  the  water  just 
beneath  the  interface  with  the  wallboard.  The  possibility  of  later  locating 
the  receiver  in  the  wallboard  is  being  considered.  The  wallboard  will  be 
separated  from  the  water  by  a  very  thin  sheet  of  plastic  so  that  the  amount 
of  water  absorbed  by  the  wallboard  can  be  controlled. 
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IV. 6  Fields  and  Currents  and  Charges  on  Obstacles  In  a  Parallel-Plate  Simulator 
at  Selected  Frequencies  and  with  Pulse  Excitation.  W.-Y.  P^^l,  R.W.P.  King 
and  T.T.  Wu,  Contract  N00014-84-K-0465;  Research  unit  11  (former  #10). 

The  parallel-plate  Hajnf-ard  BMP  simulator  has  been  investigated  during  the 
past  several  years  in  order  to  acquire  a  deep  understemding  of  its  properties 
as  a  means  to  inprove  its  quality  for  EMP  testing.  The  research  was  first 
Ccurried  out  in  the  frequency  domain  £uid  was  then  extended  to  the  time  domain. 
The  investigation  of  the  parallel-plate  simulator  under  pulse  excitation  was 
begiin  with  a  comprehensive  experimental  study  and  subsequent  theoretical  ana¬ 
lysis  of  the  sinpler  rhonibic  structure  used  as  an  EMP  simulator  excited  by  an 
inpulse  generator.  '  (These  two  papers  received  the  1963  Tremsactions  Prize 
Paper  Jlward  of  the  IEEE  Electromagnetic  Oospatibility  Society  in  October  1984.] 
An  e3q>erimental  investigation  of  the  metal-plate  EMP  simulator  under  pulse 
excitation  was  then  conducted.^  ^  The  conclusions  are:  (a)  Formulas  for  the 
electric- field  pulses  in  the  working  volume  of  the  rhombic  simulator  show  that 
these  propagate  quite  smoothly  and  uniformly.  They  are  continuously  generated 
by  current  pulses  in  the  rhombic  conductors,  (b)  The  incident  pulse  is  every¬ 
where  independent  of  the  terminating  load.  Since  in  an  EMP  test  the  incident 
pulse  alone  is  of  interest,  the  nature  of  the  termination  is  irrelevant. 

(c)  The  nature  of  the  electric-field  pulses  observed  at  successive  locations 
in  the  metal-plate  simulator  is  far  from  simple.  Their  structure  cw  be  under¬ 
stood  in  peurt  by  compeurison  with  the  corresponding  pulses  in  the  rhombic 


simulator.  Only  sliest  mismatches  between  the  pulse  generator  euid  the  coaxial 
line  and  between  this  and  the  simulator  suffice  to  tremsform  a  single  pulse 
leaving  the  generator  into  a  complex  pulse  sequence  entering  the  simulator. 

During  this  reporting  period  the  work  has  continued  with  a  study  of  the 
obstacle- simulator  interaction  under  pulse  excitation.  The  distributions  of 
surface  ciirrent  and  chuge  induced  on  a  vertical  metal  cylinder  by  am  electro¬ 
magnetic  pulse  have  been  investigated  experimentally  in  the  Harveurd  parallel- 
plate  BMP  simulator.^  Fourier  transforms  of  the  observed  pulse  shapes  cure 
taken  in  order  to  obtain  distributions  of  the  several  frequency  cooponents  on 
the  surface.  These  eu:e  then  compared  with  theoretical  distributions  calculated 
using  Kao's  method.  The  measured  data  axa  found  to  be  in  beisically  good  agree¬ 
ment  with  those  theoretical  distributions  induced  by  the  sum  of  the  incident  amd 
reflected  pulses.  Observed  differences  between  mecisured  and  theoretical 
distributions  for  certain  frequency  con^nents  of  the  transverse  current  Ig 
are  probably  due  to  coupling  between  the  cylinder  and  the  simulator  top  plate. 

A  quantitative  analysis  of  the  interaction  between  the  cylinder  and  the 
simulator  is  very  interesting  but  difficult;  it  requires  further  study. 

The  meeisured  current  pulse  sequences  induced  on  the  vertical  cylinder 
described  above  have  also  been  coit^eured,  in  a  separate  paper,  with  theoretically 
determined  pulse  shapes  of  current  induced  on  an  infinitely  long  conducting 
cylinder  by  an  E-polarized  electromagnetic  pulse.  For  an  E-poleurized  incident 
electromagnetic  pulse  with  a  nonzero  direct  component,  the  axially  induced 
current  pulse  has  an  infinite  direct  component  ajid  the  Fourier  integration  is 
divergent.  In  order  to  avoid  this  difficulty,  the  Laplace  treuisform  for  the 
variable  t  is  taken  instead  of  the  Fourier  transform.  After  the  inverse 
Laplace  transform  has  been  taken,  the  pulse  shapes  of  the  axial  induced  current 
can  be  expressed  as  the  sum  of  only  a  few  terms.  The  measured  pulse  sequences 
for  the  finite  cylinder  closely  resemble,  in  its  main  features,  those  of  the 
infinite  cylinder  but  with  a  series  of  small  reflected  pulses  added. 

The  experimental  and  theoretical  investigation  of  the  rhombic  simulator 
1  2 

under  pulse  excitation  '  indicated  that  the  electromagnetic  pulse  in  the  work¬ 
ing  volvune  of  the  simulator  is  excited  by  the  current  in  the  wires  near  the 
soiirce  rather  than  that  in  the  wires  necu:  the  center  of  the  simulator.  This 
observation  is  somewhat  different  from  that  in  the  frequency  domadn.  The  same 
question  must  be  asked  for  the  parallel-plate  simulator:  Does  the  triangular 


plate  near  the  source  also  play  the  most  importeuit  role  in  forming  the  electro¬ 
magnetic  pulse  in  the  working  volume  of  the  peurallel-plate  simulator?  To 
answer  this  question,  the  electromagnetic  pulse  has  been  evaluated  in  terms 
of  the  currents  on  the  plates  and  the  roles  of  the  currents  on  different  parts 
of  the  plates  have  been  analyzed.  A  paper  is  now  being  written  which  clarifies 
the  excitation  mechanism  of  the  parallel-plate  EMP  simulator.  It  is  fovmd  that 
if  the  length  of  the  parallel  plate  is  short  (i.e.,  close  to  the  height  h) , 
the  currents  on  the  parallel  plate  contribute  almost  nothing  to  the  electro¬ 
magnetic  pulse. 
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IV. 7  Arrays  of  Antennas  Along  Curved  Lines.  T.T.  Wu,  contract  N00014-84-K- 
0465  5md  DOE  Grant  EE-FG02-84ER40158;  Research  Unit  11  (former  #10) . 
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Over  fifty  years  ago,  Fermi  invented  his  method  of  pseudopotentials  to 
study  nuclear  scattering  as  described  by  the  three-dimensional  Schrodinger 
equation.  It  is  clearly  most  desirable  to  extend  this  method  to  radiation 
problems  2is  described  by  Maxwell's  equation.  Unfortunately,  there  is  a 
rigorous  mathematical  theorem  for  operators  in  Hilbert  space  that  such  an 
extension  is  inpossible. 

An  intermediate  problem  between  those  of  nuclear  physics  and  electro¬ 
magnetic  theory  is  described  by  the  Schrodinger  equation  in  five  dimensions. 

It  wcis  found^  that  the  rigorous  mathematical  theorem  can  be  circumvented, 
the  key  point  being  the  replacement  of  the  Hilbert  space  by  the  Pontrjagin 

space.  With  this  point  understood,  the  method  of  Fermi  pseudopotentials  has 

2 

been  generalized  to  cuitenna  2u:rays.  So  far  only  very  simple  arrays  have 
been  treated,  2md  the  other  cases  remain  to  be  worked  out  (see  p.  304  of 
ref.  2). 
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IV. 8  Solitary  Electromagnetic  Pulses  with  Slow  Decay.  T.T.  Wu  and  R.W.P.  King 
contract  N00014-84-K-0465  and  DOE  Grant  DE-FG02-84ER40158;  Resecurch  Unit 
11  (former  #10) . 

In  a  recent  paper  by  J.N.  Brittingham,^  mathematical  formulations  are 

developed  for  "new,  three-dimensional  packet-like  solutions  to  the  free-space 

Maxwell's  equations."  these  solutions  cure  further  characterized  as  "real, 

nonsingular,  continuous  functions  which  propagate  in  a  straight  line  at  light 

velocity.  They  remain  focused  for  all  time."  It  was  subsequently  shown  by 
2 

Belanger  that  a  whole  feunily  of  packet-like  solutions  of  the  homogeneous  wave 
equation  can  be  derived  that  includes  Brittingham's  form  as  a  special  case. 


IV.17 


Brittingheun  amd  Belanger  both  observed  that  pulses  of  the  type  formulated, 
while  satisfying  Maxwell's  equations,  carry  infinite  energy  so  that  they  are 
physically  unrealizable.  In  order  to  overcome  this  difficulty,  Brittingham^ 
supplemented  his  three-dimensional  wave-packet  by  preceding  and  following 
surfaces  of  discontinuity  which  served  to  keep  the  energy  associated  with  the 
packet  finite.  Work  in  this  area  at  Harvard  begem  with  the  determination^ 
that  these  bounded  pulses  of  Brittingham  no  longer  satisfy  Maxwell's  equations 
across  the  surfaces  of  discontinuity  and  aure  therefore  not  physically  meaning¬ 
ful.  Current  research  is  concerned  with  the  possibility  of  generating  electro¬ 
magnetic  pulses  which,  while  not  "focused  for  all  time"  like  Brittingham' s , 

nevertheless  propagate  in  a  manner  such  that  the  energy  associated  with  them 

-2 

decreases  with  the  radial  distance  R  much  nvore  slowly  them  as  R  .  The 
justification  for  the  belief  that  this  may  be  realizable  comes  from  recent 
advances  in  particle  physics.  In  particular,  it  is  leeurned  from  this  work 
that,  even  for  extremely  short  wavelengths,  a  logarithmic  dependence  of  such 
quantities  as  the  total  cross-section  on  the  wavelength  cam,  and  does,  occur. 

4 

In  the  specific  formulation  proposed  recently,  it  is  found  that  the  energy 

received  by  a  given  detector  must  approach  zero  as  it  is  moved  further  and 

further  away  from  the  pulse- generating  amtenna.  However,  this  approach  to 

-2 

zero  cam  be  much  slower  tham  given  by  the  fom  R  for  an  outward- traveling 
spherical  wave.  Such  cases  of  slow  decrease  are  referred  to  as  electromagnetic 
missiles.  Another  way  of  stating  the  result  is  that  the  product  of  the 
following  two  quantities  can  approach  zero  as  slowly  as  one  wishes:  (1)  the 
energy  transmitted  per  unit  aurea  of  receiver,  and  (2)  the  total  area  of  the 
receiver.  Examples  aure  constructed  when  either  one,  but  not  both,  of  these 
quamtities  remains  finite  at  infinite  distance.  Since  the  rate  of  decrease 
with  distance  of  electromagnetic  ndssiles  cam  be  very  slow  and  they  have  the 
advamtage  of  moving  with  the  velocity  of  light,  they  have  important  possible 
applications. 
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\\l.3  Advanced  Introduction  to  Electromagnetic  Theory  and  Its  Applications. 


R.W.P.  King  auid  M.  Owens,  and  S.  Prasad  Hrnchey  of  Northeaster  Univer¬ 
sity,  Contract  N00014-84-K-0465;  Research  Unit  11  (former  #10). 

The  manuscript  for  a  new  introductory  book  on  electromagnetic  theory  eind 
its  applications  has  been  con$)leted  and  is  in  press  by  Prentice-Hall.^  The 
foundations  of  electromagnetism  aire  presented  in  the  first  six  chapters  in  a 
manner  following  that  in  King's  Fundamental  Eleotromagnetio  Theory — which  is 
now  out  of  print.  Applications  begin  with  a  chapter  on  the  scattering  and 
diffraction  of  plane  waves  by  a  half-plane  and  continue  with  chapters  on 
antennas,  electric  circuit  theory,  transtnission-line  theory,  the  insulated 
antenna,  the  theory  of  metal  waveguides,  waves  along  dielectric  rods  with 
reference  to  optical  fiber  transmission,  and  electromagnetic  surface  waves 
along  boundeuries.  The  presentation  stresses  the  physical  interpretation  of 
the  mathematical  symbolism  and  the  basic  unity  of  the  theory.  The  applications 
are  arranged  with  increasing  difficulty  as  logical  continuations  of  the  theory. 
They  deal  specifically  with  proolems  of  continuing  and  current  interest. 
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V.  SIGNIFICANT  ACCOMPLISHMENT  REPORT 


V. 1  X-Ray  Reflectivity  Measurement  of  the  Structure  of  the  Surface  of  Water. 

A.  Braslau,  A.  Weiss,  J.  Als-Mielsen/  and  P.S.  Pershan,  Contracts 

N00014-84-K-0465,  NSF  DMR-82-12189 ,  NSF  DMIi-83-16979,  and  NSF  INT-83- 11841; 

Reseeurch  Unit  3. 

During  this  report  period  we  made  the  first  measurement  of  the  x-ray 
reflectivity  of  the  surface  of  water-  The  angular  dependence  of  the  reflecti¬ 
vity  is  theoretically  expected  to  follow  the  following  law : 

R(6)  =  BF (9) exp [- (Q^)a^ 

where  RF(6)  is  the  reflectivity  calculated  from  the  Fresnel  law  of  optics  on 
the  assumption  of  a  sharp  flat  surface  between  air  and  the  known  dielectric 
constant  of  water  at  1.'S'4  Angstroms ;  Q=  (47r/X)  sin  (0)  ;  and  o  is  the  mean  square 
roughness  of  the  surface .  A  simple  model  attributing  the  surface  roughness 
to  thermally  excited  capill2u:y  waves  predicts  a »  (kT/4’ry)  j  d  q/q  where  y 
is  the  surface  tension  of  water  and  the  limits  on  the  integral  depend  on  the 
spectrometer  resolution  and  the  molecular  dimensions .  The  data  agrees  with  this 
form  and  obtains  a  value  for  0  of  2.8  while  this  model  obtains  3.2.  Ihe 
dependence  on  l/y  is  confirmed  by  our  recent  measurements  on  CCl^  .  Although 
there  has  been  considerable  theoretical  research  into  the  surface  roughness  of 
simple  liquids,  this  is  the  first  measurement  on  the  molecular  length  scale. 

*This  project  was  done  in  collaboration  with  Als-Nielsen.  All  his  support 
as  well  as  the  experimental  facilities  at  the  synchrotron  ledx  in  Hamburg, 
Germeuiy  were  provided  by  the  Ris^  National  Laboratory  of  Denmark, 


Frequencies.  Q.  Hu,  M.  lemsiti,  o.  Liengme,  J.U.  Free,  euid  M.  Tinkham, 
Contracts  N00014-84-K-0465  and  N00014-83-K-0383;  Hesearch  Unit  4 
(former  #2) . 

In  last  year's  report,  we  described  our  success  in  demonstrating  that 
responsivities  near  the  quantum  limit  (one  electron  per  photon)  could  be  ob- 
tcdned  in  superconducting  tunnel  junctions  operating  in  the  far-infrared  (FIR) 
spectral  region.  The  major  difficulty  in  pursuing  this  further  in  the  direc¬ 
tion  of  practical  devices  was  that  these  experimental  junctions  were  fediricated 
from  soft  metals  (Sn  and  Pb)  which  made  them  unable  to  withstand  the  strains 
associated  with  repeated  cycling  between  room  temperature  and  helium  tenperature. 
To  address  this  problem,  we  organized  a  collaboration  with  Dr.  L.N.  Smith,  then 
of  the  Sperry  Research  Laboratories,  in  which  he  fetbricated  some  devices  using 
the  SNAP  (Selective  Niobium  Anodization  Process)  process.  Apart  from  the 
radically  different  material  and  fabrication  method,  the  configuration  of  these 

devices  was  similar  to  those  reported  earlier,  having  broadly  resonemt  half-wave 
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thin-film  antennas  and  a  junction  area  of  order  1  pm  .  Because  of  their  lower 
resistance,  larger  dielectric  constant  of  the  barrier,  and  somewhat  leurger 
physical  size,  coupling  into  these  devices  was  not  quite  as  effective  as  with 
the  home-made  Sn-Pb  junctions,  but  preliminary  experiments  showed  that  Josephson 
steps  could  be  induced.^  Moreover,  because  of  the  lower  resistance,  noise 
rounding  effects  were  minimal  compeured  to  those  in  the  higher-resistamce  Sn-Pb 
junctions.  The  shauT)er  resulting  features  should  be  superior  for  detection 
purposes,  which  depend  on  these  extreme  nonlineeurities . 

In  the  process  of  exploring  the  properties  of  these  junctions,  we  dis¬ 
covered  that  merely  reducing  the  FIR  laser  frequency  from  604  GHz  to  419  GHz 

(or  to  245  GHz)  resulted  in  observation  of  I-V  curves  of  an  entirely  different 
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character  than  those  seen  before,  emd  so  thoroughly  understood;  we  attribute 
these  radical  changes  to  the  onset  of  "chaos"  in  the  driven  motion  of  the 
junction.  The  hallmaurlcs  observed  aure:  extremely  high  levels  of  low-frequency 
noise,  subharmonic  step  structures  even  when  adjacent  integer  steps  are  missing, 
emd  meandering  voltage  levels  in  the  I-V  curve  on  otherwise  well-defined  and 
flat  steps.  Moreover,  these  effects  are  most  prominent  when  the  FIR  drive 
frequency  is  close  to  the  Josephson  plasma  frequency  (estimated  at  420  G&z  in 
these  junctions) ,  as  is  expected  from  computer  sithulation  studies  of  chaos. 


The  low-frequency  noise  levels  in  the  chaotic  regions  can  be  as  large  as 

the  equivalent  of  lO^^K  thermal  noise  at  100  Hz-  Further  measurements  of  the 
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power  spectrum  show  that  it  falls  off  as  1/f  at  higher  frequencies,  where 
a  varies  from  1  to  4,  and  may  depend  on  the  levels  of  dc  bias  and  laser  drive 
power.  On  the  basis  of  our  simulations,  we  suggest  that  although  the  noise 
levels  may  be  qualitatively  similar  in  different  operating  points,  the  fvinda- 
mental  origins  may  involve  different  chaotic  mechanisms,  such  as  period  doub¬ 
ling  cascades,  intrinsic  intermittency ,  and  noise- induced  intermittency . 

One  region  we  have  studied  particularly  closely  is  that  at  the  first 
onset  of  finite  dc  voltage,  in  which  the  system  appears  to  be  intermittently 
shifting  between  a  zero-step  and  2/3  step  regime  with  a  "1/f"  spectrum.  Our 
conputer  simulations  (on  a  VAX)  show  that  the  "basins  of  attraction"  for  these 
two  steps  are  extremely  conplex  in  form,  separated  by  a  boundary  having  a 
"fractal"  dimension  which  depends  on  dc  bias  level  and  can  grow  almost  to  2  (at 
which  point  the  boundary  would  fill  the  entire  phase  pl^me} .  The  experimental 
noise  level  shows  a  dependence  on  bias  current  that  closely  parallels  this 
confuted  fractal  dimension,  qualitatively  confirming  the  notion  of  Grebogi, 
et  al.^  that  these  fractal  dimensions  imply  a  sensitivity  to  initial  conditions 
(emd  hence  noise)  which  diverges  as  the  dimension  approaches  2.  Effectively, 
shot  noise  due  to  the  discreteness  of  the  electronic  cheurge  is  an^lified 
enormously  to  generate  immense  1/f-type  noise.  This  work  is  the  first  of  which 
we  are  awcure  in  whidi  data  on  a  real  physical  system  has  confirmed  these  isplicit 
consequences  of  the  fractal  basin  boundcury.  The  practical  significance  of  this 
work  is  that  the  numy  routes  to  chaos  must  be  understood  so  that  this  source  of 
extraordinary  noise  levels  can  be  avoided  in  the  design  of  actual  devices. 
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Unit  11  (former  #10). 


An  importwt  application  of  electromagnetic  waves  is  in  remote  sensing, 
notably  in  the  search  for  underground  regions  containing  oil  or  ore.  A  con^ 
venient  and  useful  procedure  makes  use  of  a  fixed  insulated  transmitting  dipole 
and  a  movable  receiving  antenna  both  on  the  surface  of  the  earth.  How  is  the 
field  measured  by  the  receiver  related  to  a  buried  object?  A  quantitative 
cuiswer  to  this  question  is  essential  to  an  understanding  and  eiqploitation  of 
the  method  and  the  correct  interpretation  of  measured  data.  This  has  now  been 
provided  in  the  form  of  a  con^lete  theoretical  analysis  emd  a  significeuit  set 
of  preliminary  measurements  in  a  model  t£uik. 

The  first  step  in  the  analysis  is  the  determination  of  the  electromagnetic 

field  incident  on  a  buried  object.  This  is  not  the  field  that  travels  directly 

through  the  esirth  from  the  trcuismitter  to  the  scattering  obstacle  or  boundary 

but  a  lateral  wave  that  proceeds  along  the  surface  in  the  air  to  a  point  above 

the  obstacle  and  then  vertically  dcwn  to  it.  Recent  reseeurches  have  provided 
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simple  auid  accurate  formulas  for  the  six  components  of  this  field.  They 

have  been  determined  specifically  for  a  horizontal  dipole  over  the  earth. 

Signif iceuitly ,  the  largest  con^onents  of  this  field  are  the  radial  electric 

(Ej^p)  euid  tr6msverse  magnetic  direction  along  the  axis  of  the 

trauismitting  antenna.  In  order  to  obted.n  a  single  and  accurate  solution  that 

could  be  used  to  verify  an  experiment  in  a  model  tank,  the  first  choice  of 

buried  c^ject  to  be  located  was  an  insulated  cylindrical  conductor  in  salt 

water.  Oonplete  calculations  of  the  incident,  scattered,  and  total  fields 

over  a  wide  area  over  the  surface  of  the  earth  above  the  submerged  cylinder 

with  its  depth  as  parameter  show  a  clear,  readily  interpreted  interference 
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pattern  directly  over  the  cylinder.  The  location,  shape  and  orientatxon  of 
the  cylinder  sure  clearly  determined  as  sharp  rises  or  dips  in  the  total  field. 
This  interference  pattern  consists  of  sharp  increases  in  the  field  over  the 
cylinder  when  its  depth  is  such  that  the  scattered  emd  incident  fields  are  in 
phase  at  the  receiver,  or  equally  sharp  decreases  when  the  incident  emd 
scattered  fields  are  out-of-phase. 


NO  measurements  of  the  scattered  field  of  a  submerged  metal  cylinder 
have  yet  been  made.  However,  measurements  have  been  carried  out  by  Bansal 


of  the  corresponding  patterns  generated  by  a  submerged  circular  metal  disk. 

These  indicate  the  same  type  of  interference  above  the  edges  of  the  disk.  To 
complete  the  study,  measurements  should  be  made  for  the  suibmerged  cylinder  and 
the  theory  developed  for  the  submerged  disk.  Besults  to  date  fully  confirm 
and  illuminate  the  method  and  describe  the  conditions  for  its  optimum  applica¬ 
tion.  An  even  broader  scc^e  could  be  achieved  with  theoretical  studies  and 
sets  of  measurements  for  a  number  of  different  scattering  regions  with  proper¬ 
ties  other  than  those  of  conductors. 
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